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rNTRODUCTION 


The  present  investigation  was  carried  out  by 
Dr«  K»  0.  Prledrlchs  of  the  Institute  for  Mathematics  and 
Mechanics  at  Hew  York  University  under  Contract  N0a(s)-7370, 
Bureau  of  Ordnance,  Navy  Department*  It  continues  the  work 
done  for  the  Applied  Mathematics  Panel  of  the  National  Defense 
Besearch  Committee  In  particular  the  Supplement,  "Remarks  on 
the  Mathematical  Theory  of  Detonation  and  Deflagration  Waves 
Ijx  Gases,"  to  the  Manual  on  "Supersonic  Plow  and  Shock  Waves," 

Thanks  are  due  to  Dr.  Robert  Fli^elsteln  and  Dr*  James 
Savage  for  many  critical  comments  and  sxiggestlons. 


R,  Courant 
Director 
laatltute  for  Mathematics  and  Mechanics 


On  the  Mathematical  Theory  of  Deflagrations 
and  Detonations 

While  the  propagation  of  a  shock  wave  Is  coiDpletely  deter- 
Bdned  by  the  conservation  laws,  the  boundary  conditions  of  the 
problem^  and  the  additional  condition  that  the  entropy  Increase 
in  the  process,  the  same  is  not  true  for  the  propagation  of  a 
detonation  wave  and  of  the  flame  front  in  an  ordinary  combustion 
process •   More  conditions  must  be  added  to  the  conservation 
laws  in  order  to  provide  sufficient  data  for  the  unique  deter- 
mination of  the  propagation  process.  For  detonations  this 
necessity  was  recognized  by  Chapman  and  Jouguet  when  they  intro- 
duced their  famous  hypothesis*   For  combustion  processes  this 
necessity  was  more  or  less  tacitly  assumed  by  Jouguet  and  others 
when  they  attacked  the  calculation  of  the  flame  speed  by  taking 
heat  conduction  into  account  without  even  trying  to  determine 
the  flame  speed  from  the  conservation  laws  and  boundary  condi- 
tions alone* 

It  is  natural  to  expect  that  the  needed  additional  condi- 
tions could  be  derived  from  an  investigation  of  the  Internal 
mechanism  of  the  combustion  or  detonation  process.   Thus  v«  Heunaxm 
[l]  has  arrived  at  a  justification  of  the  Chapman- Jouguet  hypoth- 
esis for  detonations  by  taking  into  account  that  the  chemical 
reaction  takes  plabe  over  a  zone  of  finite  width;  his  arguments 
are  based  on  the  assus^tion  that  a  detonation  is  initiated  by  a 
shock.  For  combustion  processes,  which  do  not  involve  shocks, 
no  unique  determination  can  be  achieved  without  taking  heat  oon- 
duction  into  account. 

It  is  the  intention  of  the  present  paper  to  offer  a  unified 
and  more  complete  discussion  of  the  question  of  determinacy  for 


detonations  and  deflagrations,'   In  order  to  be  able  to  point 
out  the  contrast  between  these  two  kinds  of  processes  we  shall 
treat  both  of  them  on  the  basis  of  the  same  assumptions:  We 
shall  take  viscosity  and  heat  conduction  into  accoimt  and  assume 
a  finite  rate  of  chemical  reaction.   (Accordingly,  we  shall  not 
postTilate  that  a  detonation  process  begins  with  a  shock.)  Prom 
the  discussion  of  the  internal  mechanism  of  the  detonation  and 
deflagration  process  on  this  basis  we  shall  obtain  the  desired 
additional  conditions  which  make  unique  determination  of  the 
whole  process  possible  by  excluding  certain  detonation  or  defla- 
gration processes  which  would  be  compatible  with  the  conservation 
laws.  In  particular  we  shall  find  as  a  result  that  a  detonation 
begins  with  a  shock  smd  that  the  Chapman-Jouguet  hypothesis 
furnishes  the  correct  additional  condition  provided  that  for  a 
given  value  of  the  reaction  rate  the  viscosity  and  the  heat  con- 
ductivity are  sufficiently  small.  If,  on  the  other  hand,  the 
reaction  rate  is  very  high,  for  given  viscosity  and  heat  conduc- 
tivity, the  detonation  no  longer  begins  with  a  shock;  and  if  the 
reaction  rate  is  excessively  high,  the  Chapman-Jouguet  hypothesis 
is  no  longer  correct;  (see  Appendix  II),  The  additional  condi- 
tion for  deflagrations  is  such  that  it  prescribes  the  flame  speed 
dependent  on  heat  conduction  and  reaction  rate. 

To  fix  the  ideas  we  have  considered  the  problem  of  determin- 
ing the  flow  under  the  following  circumstances:   In  an  Infinite 
tube  a  piston  moves  in  a  prescribed  manner  beginning  at  an  Initial 
time  at  an  initial  cross-section.  At  the  same  time  a  reaction 
process  begins  at  the  piston  and  travels  into  the  quiet  unburnt  gas, 

1)  We  propose  to  use  the  term  deflagration  for  those  combustion 
processes  which  take  place  in  a  very  narrow  zone  of  constant  width 
and  which  therefore  in  good  approximation  can  be  described  by  a 
discontinuity.  For  detonation  and  deflagration  processes  we  shall 
employ  the  common  name  "reaction  process". 


This  problem  Includes  the  Important  case  that  the  tube  is 
closed  at  the  initial  cross-section;  for,  this  case  restilts 
when  the  piston  remains  at  its  original  place.  The  case  of 
an  open  end  results  when  the  piston  is  withdrawn  with  suffi- 
ciently large  speed.  Also  the  case  is  included  that  in  a 
doubly  infinite  tube  containing  combustible  and  non-combustible 
gas,  not  separated  by  a  piston,  a  reaction  begins  at  the  inter- 
face. One  needs  only  force  the  piston  to  move  in  the  same 
manner  in  which  the  initial  cross-section  would  move  anyhow  if 
there  were  no  piston. 

Before  describing  the  results  in  greater  detail  we  formu- 
late the  basic  assximption  underlying  our  investigations  viz, 
that  the  reaction  process  in  question  may  approximately  be  con- 
sidered as  a  sharp  discontinuity.  More  precisely  the  assumption 
is  that  the  "reaction  zone,"  across  which  chemical  composition, 
pressure,  temperature,  and  velocity  change  is  very  narrow  and 
of  nearly  constant  width;  (see  Appendix  I),  By  this  we  mean, 
firstly,  that  the  space  rates  of  change  of  the  pertinent  quanti- 
ties over  a  section  in  the  field  of  flow  outside  of  the  reaction 
zone  are  negligibly  small  when  compared  to  the  rates  of  change 
of  the  same  quantities  inside  the  reaction  zone  and,  secondly, 
that  the  time  rate  of  change  of  the  width  of  the  reaction  zone 
is  small  when  compared  with  the  average  speed  with  which  the 
gases  cross  the  reaction  zone.  This  assumption  appears  to  be 
Justified  only  if  the  coefficients  of  viscosity  and  heat  con- 
duction are  sufficiently  small  and  the  rate  of  reaction  is  suffi- 
ciently large,   (We  leave  aside  the  question  whether  or  not  under 
these  circumstances  the  assumption  is  always  Justified.)  The 
pertinent  quantities  on  both  sides  of  the  "reaction  front"  are 
then  connected  by  the  same  well-known  laws  of  conservation  of 
mass,  momentum,  and  energy  that  hold  for  the  quantities  at  both 
sides  of  a  discontinuity  surface. 

Next  we  give  a  brief  accoiont  of  the  indeterminacies  that 
one  encounters  when  one  tries  to  determine  a  flow  involving  a 


reaction  discontinuity  solely  by  using  the  conservation  laws 
and  the  boundary  conditions.  To  this  end  It  Is  necessary  to 
distinguish  various  types  of  reaction  processes.   Among  the 
detonations  there  Is,  as  Is  well  known,  a  particular  one,  the 
" Chapman- Jouguet"  detonation,  which  Is  singled  out  from  others 
by  the  property  that  the  flow  of  the  burnt  gas  Is  sonic  when 
observed  from  the  reaction  front.  We  have  termed  "strong"  or 
"weak"  a  detonation  If  it  Involves  a  pressure  rise  greater  or 
less  than  for  a  Chapman- Jouguet  detonation.    Similarly,  we 
have  termed  "strong"  or  "weak"  a  deflagration  If  It  Involves 
a  pressure  drop  greater  or  less  than  for  a  Chapman- Jouguet  de- 
flagration, which  again  Is  characterized  by  the  condition  that 
the  flow  burnt  gas  Is  sonic  when  observed  from  "the  reaction 
front.  A  constant  volume  detonation  Is  the  limiting  case  of 
a  weak  one,  producing  the  least  pressure  rise  (and  the  least 
tenperature  rise)  among  all  detonations  for  the  same  explosive, 
a  constant  pressure  deflagration  Is  also  the  limiting  case  of 
a  weak  one,  producing  the  least  drop  In  density  (and  the  greatest 
temperature  rise)  among  all  deflagrations. 

We  now  consider  the  flow  of  gas  In  a  half -finite  tube 
resulting,  as  Indicated  before,  when  a  reaction  front  starts  to 
move  from  the  finite  end  of  the  tube  into  the  unburnt  gas  under 
the  influence  of  a  piston  which  moves  In  a  prescribed  manner. 
We  then  ask  for  the  gas  motions  which  are  compatible  with  the 
conservation  laws  and  the  piston  motion.   Mathematically  speaking, 
we  ask  for  the  solutions  of  the  flow  differential  equations  com- 
patible with  the  transition  conditions  at  the  discontinuity  front 
and  with  the  boundary  condition,  which  expresses  that  the  gas 
adjacent  to  the  piston  has  the  same  velocity  as  the  piston.   This 
problem  will  be  referred  to  as  the  "external  flow  problem."   The 

1)  For  the  following  see  the  Supplement  [3]  to  the  Manual  on 
Supersonic  Plow  and  Shock  Waves  [2]  . 


answer,  explained  In  detail  In  an  earlier  report  [3]  Is  this: 
Suppose  the  reaction  process  Is  a  detonation;   If  then  the  piston 
moves  In  the  same  direction  as  the  reaction  front,  and  If  the 
velocity  Up  of  the  piston  exceeds  the  gas  velocity  Ujj  which 
would  be  produced  by  a  Chapman- Jouguet  detonation,  then  there 
Is  just  one  flow  Involving  a  strong  detonation  In  which  the 
burnt  gas  has  the  prescribed  piston  velocity.  If,  however,  the 
piston  velocity  Up  Is  less  than  Ujj,  adjustment  of  the  velocity 
of  the  burnt  gas  to  the  piston  velocity  can  always  be  achieved 
by  a  Chapman- Jouguet  detonation  followed  by  an  appropriate  rare- 
faction wave,  but  It  can  also  be  achieved  by  a  set  of  weak  detona- 
tlons  followed  either  by  a  shock  or  a  rarefaction  wave*   Thus,  In 
case  Up  <  Ujj,  the  solution  of  the  mathematical  problem  Is  not 
unique;  there  Is  a  one -parametric  set  of  solutions. 

For  deflagration  processes  the  degree  of  Indeterminacy  1b 

still  higher.  To  any  piston  velocity  Unone  has  still  the  choice 

z) 

of  a  flame  velocity  arbitrary  within  certain  limits  '  and  can 

achieve  adjustment  of  the  gas  velocity  to  the  piston  velocity 
by  sending  ahead  of  the  flame  a  shock  of  appropriate  strength. 
Again,  there  Is  a  one -parametric  set  of  solutions  as  long  as  the 
deflagration  remains  weak.  If  the  resulting  deflagrations  be- 
come strong  ones,  the  possibilities  of  adjustment  become  still 
greater  and  the  set  of  solutions  Is  two-parametric. 

To  express  these  determlnacy  statements  In  a  slinple  form 
we  Introduce  as  "degree  of  under-determlnacy"  of  the  external 
flow  problem  the  number  t  of  conditions  that  mast  be  Inposed  on 

1)  It  should  be  emphasized  that  the  theory  as  considered  In  this 
report,  based  on  the  assuirptlon  that  the  reaction  front  be  a 
sharp  discontinuity,  does  not  offer  any  possibility  of  predicting 
whether  a  detonation  or  a  deflagration  will  occur  in  a  given 
situation  (except  that  under  certain  circumstances  deflagration 
flow  is  not  possible). 

2)  The  limiting  case  is  the  one  in  which  the  flame  together  with 
the  pre-conpression  shock  are  just  equivalent  to  a  detonation. 


the  data  of  the  flow  problem  In  order  to  make  the  solution 
unique.   Summarizing  we  then  have: 

Strong  detonations      t  =  0  , 

Weak  detonations        t  =  1  , 

Weak  deflagrations      t  =  1  , 

Strong  deflagrations    t  =  2  . 
Chapman- Jouguet  detonations  and  deflagrations  are  here  classed 
with  strong  detonations  or  weak  deflagrations  respectively. 

These  peculiar  under-determinacies  are  a  consequence  of 
Jouguet' 3  ingjortant  rule  (cf.[3])j  concerning  the  properties  of 
the  gas  flow  observed  from  the  reaction  front: 
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Figure   1 


The   flow  ahead  of  a  detonation  front  la  supersonic ,. 
and  subsonic  ahead  of  a  deflagration  front*   The  flow  Is  sub- 
sonic behind  a  strong  detonation  and  a  weak  deflagration  front. 
supersonic  behind  a  weak  detonation  and  a  strong  deflagration 
front*   The  flow  is  here  always  understood  relative  to  the 
reaction  front* 

We  now  proceed  to  discuss  the  chief  aim  of  the  present 
paper,  namely  to  decide  which  of  the  flow  processes,  still  per- 
mitted by  the  conservation  laws,  are  excluded  through  the  action 
of  viscosity,  heat  conduction,  and  chemical  reaction.   To  this 
end  we  shall  take  into  account  that  the  reaction  zone  has  a 
finite  extension.  We  then  shall  set  up  the  differential  equa- 
tions governing  the  transition  across  such  a  reaction  zone  and 
investigate  under  which  circumstances  these  differential  equa- 
tions possess  solutions  satisfying  the  boundary  conditions  im- 
posed at  the  two  ends  of  the  reaction  zone.   These  boundary 
conditions  consist  in  prescribing  the  chemical  composition, 
pressure,  ten5)erature,  and  velocity  of  the  gases  at  both  ends 
of  the  reaction  zone  in  such  a  way  that  the  laws  of  conservation 
of  mass,  momentum,  and  energy  are  satisfied  by  these  quantities. 
The  differential  equations  in  the  interior  of  the  reaction  zone 
express  the  same  conservation  laws,  but  take  into  account  chemi- 
cal reaction,  viscosity,  and  heat  conduction. 

Investigations  of  this  kind  for  shocks  not  Involving  a 
chemical  reaction  have  been  made  in  detail  by  various  authors  , 
The  result  was  that  a  transition  between  the  given  quantities 
at  both  sides  of  the  zone  is  always  possible  provided  that  the 
direction  of  the  flow  corresponds  to  increasing  entropy.   Our 
investigation  will  show  that  the  same  is  true  for  strong  detona- 
tions and  for  Chapman- Jougue t  detonations.   For  oth»p  modes  of 
reaction,  however,  the  situation  is  completely  different.  Here 
are  the  results  of  our  analysis: 

1)  See  e.  g.  Becker  [4]  and  Weyl  [5]  • 


Unless  the  rate  of  reaction  is  excessively  high,  weak 
detonations,  though  compatible  with  the  conservation  laws, 
are  impossible.  The  condition  that  the  detonation  be  strong 
or  of  the  Chapman- Jouguet  type  is,  therefore,  the  desired  addi- 
tional condition  mentioned  in  the  beginning.   Consequently,  a 
flow  involving  a  detonation  is  \iniquely  determined.  For,  if 
the  piston  velocity  is  high,  Up  >  vu,    there  is  a  \mique  solution 
involving  a  strong  detonation,  as  mentioned  before;  for  lesser 
piston  velocity,  however,  Up  5  vu,    there  is  now  only  one  possible 
flow  left  in  which  the  velocity  of  the  burned  gas  equals  the 
piston  velocity,  and  that  is  the  flow  involving  a  Chapman- Jouguet 
detonation.   In  particoilar  we  see  that  for  a  tube  with  a  closed 
end.  Up  =  0,  or  open  end,  Up<<  0,  a  Chapman- Jouguet  detonation 
is  the  only  possible  one.  Thus  the  occ\arrenc^  of  this  particular 
detonation  is  here  deduced  and  no  additional  hypothesis  is  re- 
quired. 

If,  however,  the  reaction  rate  is  excessively  high,  the 
analysis  yields  the  result  that  the  Chapman-Jouguet  detonation 
is  impossible.   Instead,  a  particular  weak  detonation  is  possible 
which  travels  with  a  well-determined  velocity  (depending  on 
pressure  and  temperature  in  the  lonburnt  gas,  reaction  rate,  viscos- 
ity, and  heat  conductivity).   If  the  piston  velocity  is  large 
enough,  adjustment  still  requires  a  strong  detonation.   For  lesser 
piston  velocities,  for  example  for  closed  or  open  ends,  adjust- 
ment is  effected  by  the  particular  weak  detonation  followed  by  an 
appropriate  shock  or  an  appropriate  rarefaction  wave. 

As  to  deflagrations  the  results  of  the  analysis  are  that 
strong  deflagrations  are  impossible  altogether.  '  A  weak  deflagra- 
tion is  possible  only  v;ith  a  well-determined  speed.  We  distinguish 
the  "flame  velocity,"  i.e.  the  velocity  of  the  reaction  front 
relative  to  the  tube,  from  the  "burning  speed,"  i.e.  the  speed  of 

1)  This  result  could  also  be  established  by  v.  Neumann's  argument 
Ignoring  viscosity  and  heat  conduction  and  taking  only  the  finite 
rate  of  chemical  reaction  into  account. 


the  reaction  front  relative  to  the  unburnt  gas  ahead  of  it. 
While  the  flame  velocity  depends  on  the  boundary  conditions 
of  the  problem  as  a  whole,  the  burning  speed  depends  only 
on  pressure  and  ten^erature  in  the  unburnt  gas,  and  also  on 
reaction  rate,  heat  conductivity  and  viscosity*   That  the 
burning  speed  has  this  particular  value  is  the  desired  addi- 
tional condition  for  deflagrations.  As  stated  before,  for  a 
burning  speed  arbitrary  within  certain  limits,  a  deflagration 
flow  can  be  found  which  is  adapted  to  the  piston  motion.   A 
further  limitation  is  iiiposed  by  excluding  strong  deflagrations. 
Thus  we  see:  Within  certain  limits  for  the  data  of  the  problem 
there  exists  a  uniquely  determined  flow  involving  a  deflagration 
and  adapted  to  the  piston  motion. 

We  call  attention  to  a  number  of  detailed  Investigations 
of  the  transition  process  in  reaction  zones.  Detonation  transi- 
tions have  been  determined  in  detail  by  Eyring  and  his  collabora- 
tors .  Deflagration  flame  speeds  have  been  calculated  by  Jouguet 

2) 
and  others  ,  on  the  basis  of  various  assun5)tlons  about  the  de- 
tails of  the  transition  process.   The  aim  of  the  present  report 
Is  different:   it  is  not  intended  to  give  new  methods  for  cal- 
culating such  transition  processes.  The  intention  is  solely 
to  decide  the  question  of  determinacy  by  investigating  such  transi- 
tion processes  systematically  on  the  basis  of  the  existing  theory. 


1)  See  [6]  for  an  excessively  high  reaction  rate,  [7]  for  actual 
reaction  rates  •  the  latter  paper  contains  a  great  number  of 
detailed  investigations  concerning  detonations  primarily  of 
solid  explosives. 


2)  Cf.  Lewis  and  v.  Elbe  [8],  Jost  [9],  Semenov  [10],  further 
[11]  and  [12].   the  latter  report  gives  a  survey  of  earlier  work. 
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1,  The  Internal  mechanism  of  the  reaction  process. 

We  assume  the  process  to  be  strictly  one-dimensional  ' 
and  observe  the  process  from  a  frame  moving  with  Instantaneous 
flame  or  detonation  speed.  We  assume  that  the  flow  so  observed 

is  steady  at  the  time  considered  in  the  neighborhood  of  the 

2) 
reaction  zone   ,  This  is  in  agreement  with  our  basic  assumption 

(p.  3)  that  the  reaction  zone  is  very  narrow  and  of  nearly  con- 
stant width.   Consequently,  all  quantities  depend  only  on  an 
abscissa  x,  and  not  on  the  time.   At  each  place  x,  there  is  a 
mixture  of  bvirnt  and  unburnt  gas;  we  denote  by  £,  the  fraction 
of  mass  of  burnt  gas  in  the  mixture.  We  denote  pressxire  and 
specific  volume  by  p  and  t  and  introduce  the  "reduced  temperature" 
0  =  pT  ,  a  quantity  which  has  the  dimension  of  velocity  squared. 

We  ass\ame  bxirnt  and  unburnt  gas  to  be  ideal.  Accordingly,  9  is 

3) 
proportional  to  the  temperature.  The  internal  energy  e  per  unit 

mass  of  the  burnt  and  of  the  unburnt  gas  is  assumed  to  be  a 

function  of  ©  only,  (actually  e  depends  also  somewhat  on  the 

3) 
pressvire  p,  see  also  footnote   ;  denoting  by  g  the  energy  of 

formation  per  unit  mass  (at  absolute  zero  temperature),  we  intro- 
duce the  total  energy  per  unit  mass 

E  =  e  +  g  , 
Burnt  and  \xnburnt  gas  are  distinguished  by  the  superscripts  (1) 
and  (0).  The  total  energy  per  unit  mass  of  the  mixture  is  then 

1)  This  assvimption  is  no  serious  restriction  for  the  discussion 
of  the  Internal  mechanism  since  the  flow  in  the  neighborhood  of 
a  point  of  a  reaction  front  may  be  considered  one-dimensional 
to  the  same  degree  of  accuracy  as  the  assviraptlon  is  valid  that 
the  reaction  front  is  a  sharp  discontinuity. 

The  assumption  of  the  one-dimensional  character  of  the  flow  is, 
however,  a  serious  restriction  for  the  external  flow  problem  since 
this  assumption  is  known  to  be  never  quite  satisfied  for  combustion 
waves, 

2)  The  assumption  of  "local  steadiness"  does  by  no  means  imply  that 
we  consider  only  reaction  processes  that  proceed  with  constant 
velocity  into  the  unburnt  gas.   See  the  discussion  in  the  Appendix. 

3)  The  absolute  temper atiire  is  given  by  R©/m  where  R  is  the  gas 
constant  and  M  the  molecular  weight.  We  disregard  the  dependence 
of  the  molecular  weight  on  the  mixture  ratio  £  . 


n 


On  the  energy  fimctions  E(e)  we  require  that 

II    <  ° 

or 

This  condition  implies  that  the  "liberated  energy" 

E^°^0)  -E^Ho) 

is  positive,   (It  is  convenient  to  make  this  assiimption  although 
most  of  oxor  conclusions  hold  without  it.)   Our  requirement  is 
satisfied  if  burnt  and  \mburnt  gas  are  "polytropic";   i,  e,  if 
the  energies  are  given  by 

with  constant  "^  ,  ^^,  g  ,  g^,  provided  that  the  temperatxire  is 
below  a  certain  limit  ^    • 

By  V  we  denote  the  velocity  of  the  steady  gas  flow  and  by 
m  =  T"-^v  the  "mass  flux"  of  mixture  through  a  unit  cross-section 
per  unit  time.  By  S  we  denote  the  mass  of  burnt  gas  created  per 

'This  limit  corresponds  to  ©  ^ — 2 ^ (g^-  g^).  If  one 

«o  "  ''I 
assximes  ^     =  1,4,  ^- -^1,2,  the  molecular  weight  M  =  M^-^SO,  and 
the  "liberated  energy"  g  -  g-,'^»7  kcal/gr  ,  then  the  limit  is 

about  3200°K,  This  case  is,  however,  \xnrealistic  since  for  such 
temperatxires  the  value  of  V  for  combustible  gases  will  hardly 
ever  be  as  high  as  1,4, 
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unit  mass  of  unburnt  gas  per  xinlt  time.  We  assume  that  the 
"reaction  rate,"  S,  depends  on  9  and  p,  (cf.  footnote  l)  on  p. 11), 
and  that  S  vanishes  below  a  certain  "safety  temperature": 

S  =  0  for  9  <  9g.^^ 

The  transition  between  the  two  states  on  both  sides  of  the 

reaction  front  is  effected  by  the  action  of  viscosity  and  heat 

2) 
conduction   .   We  introduce  corresponding  coefficients  jjl  and  X 

such  that 


-^^ 


ana--^^ 


are  the  viscous  pressure  and  heat  per  unit  mass  conducted  through 
a  unit  cross-section  per  unit  time  •   These  coefficients  depend 
on  p,  T  ,  and  e;  but  we  need  pay  only  little  attention  to  this 
dependence  for  our  general  discussion. 

1)  The  latter  assunption  is  made  only  to  achieve  mathematical 
simplicity.   In  recent  papers  [7],  [11],  [12],  the  reaction  rate 
is  assumed  to  be  of  the  form 

-A/e 
S  =  Sooe 

where  A,  proportional  to  the  activation  energy,  is  so  large  that 
S  is  negligible  when  9  assumes  values  corresponding  to  a  tempera- 
ture of  3000K.   The  reduced  safety  temperature  9^  has  no  precise 
significance  (as  the  ignition  tenperature  had  in  the  older  litera- 
ture); it  is  siii5)ly  a  value  below  which  S  can  be  set  equal  to  zero 
for  all  practical  purposes.)   The  maximal  reaction  rate,  S«,,  may 
depend  on  P. 

2)  We  ignore  diffusion  and  radiation.  Diffusion  should  not  be 
neglected  in  the  actual  calculation  of  flame  speeds  according  to 
[12].   We  felt  that  for  the  sake  of  sinpllcity  we  could  disregard 
diffusion  since  the  additional  terms  due  to  it  would  not  seem  to 
entail  any  essential  change  in  the  general  results. 

One  might  object  to  using  the  notions  viscosity,  heat  conduction, 
and  diffusion  if  the  width  of  the  transition  zone  is  extremely  small. 
It  seems  likely,  though,  that  nevertheless  our  results  remain  cor- 
rect in  qualitative  respects,  in  particular  as  far  as  determlnacy 
is  concerned. 

3)  The  customary  coefficients  of  viscosity  and  heat  conduction  are, 
in  our  notation,  3/4yiA  ana  A/R»  R  being  the  gas  constant  such  that  R9 
is  the  tenperature. 
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We  now  formulate  the  lav^rs  governing  the  process.  The 
continuity  equation  simply  assumes  the  form 

(0)  m  =  const. 

The  law  of  conservation  of  momentum  is 

(1)'  "/"S  +  P  +  mv=  P  =  const. 

Conservation  of  energy-^'  is  expressed  by 

(2)'  -   X||+Bip)(0)    +|v2]+v[p-^f]   =mQ=  const. 

The  balance  between  burnt  and  unburnt  gas  is  given  by 

(3)'  -  ^  H  •"  (1-^)  S(0)  =  0, 

2) 
assuming  a  first  order  reaction  from  a  unimolecular  mechanism  '  , 

The  problem  is  to  investigate  possible  solutions  of  these 
differential  equations  if  the  values  of  the  quantities  v,p,  t, 
and  <£  are  given  at  the  end  points  of  the  reaction  zone.  We 
modify  this  problem  by  prescribing  the  same  values  of  these  quanti- 
ties at  X  =  00  and  x  =  -  oo  .   That  it  is  justified  v/ith  good  approxi- 
mation to  substitute  the  modified  problem  for  the  original  one 
follows  from  our  basic  assumption  (p.  5  )  that  the  v/idth  of 

1)  If  we  were  to  consider  diffusion  we  would  introduce  a  coefficient 
S  such  that  -  %   de/dx  is  the  fraction  of  mass  of  burnt  gas  diffus- 
ing through  a  cross -section  per  unit  time.   Then  we  'would  have  to 
add  the  term  Td(  Sd5/dx)dx  to  equation  (5)'  and  -  Sd[E(£)(9)+l  v2]/dx 

to  equation  (2)'  in  order  to  express  the  diffusion  of  enerQr  .  The 
resulting  modified  equation  (2)'  would  differ  somewhat  from  those 
in  the  literature  (see  [12]  where  only  the  diffusion  of  the  energy 
of  formation  is  taken  into  account. 

2)  If  a  different  mechanism  of  reaction  were  assumed  leading  to 
terms  (1  -  &)2   s  or  (1  -  £)(1  +  ne)S,  cf.  [11]  and  [12],  no  change 
in  the  general  conclusions  would  result. 
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the  reaction  zone  is  very  narrow.  More  specifically,  the 
assumption  was  that  the  rate  of  change  of  the  pertinent  quanti- 
ties outside  the  reaction  zone  is  negligibly  snail  as  compared 
to  the  rates  of  change  of  these  quantities  inside  the  reaction 
zone.   Consequently,  these  quantities  appear  to  be  nearly  constant 
at  the  ends  of  the  reaction  zone  over  a  region  whose  extension  is 
large  compared  with  the  width  of  the  reaction  zone.   It  is  then 
natural  to  assume  that  the  process  inside  the  reaction  zone  can 
very  well  be  approximated  by  a  process  that  extends  over  the  whole 
field  from  x=-a>tox=+oo  and  in  which  the  pertinent  quanti- 
ties assiome  at  ±  oo  those  values  that  are  prescribed  for  the  pro- 
per process  at  the  ends  of  the  finite  reaction  zone.^^ 

Accordingly  we  ask  for  solutions  of  the  equations  (1)  ,(2)  , 
and  (3)'  which  are  defined  for  -  oo  <  x  <  oo  and  which  approach 
finite  limit  values  (with  Tf^O)  asx->iao;  then  the  deriva- 
tives approach  zero  as  seen  from  (1)',  (2)',  and  (3)  .  Solutions 
which  behave  in  that  way  at  x  =  oo  or  x  =  -  as  will  be  called 
"regular"  there.  The  limit  values  for  x  -*►  -  oo  are  denoted  by  p^, 
Tq,  Qq,  and  Sq*  those  for  x  -s>  ooby  p^,  T^,  ©^,  and  £^.  The 
boundary  conditions  then  consist  in  prescribing  these  values;  in 
partlculsr  we  prescribe 

expressing  that  the  gas  consists  of  unburnt  gas  at  x  =  -  co  and 
of  completely  burnt  gas  at  x  =  +  oo  .  Prom  relation  (5)'  we  then 

1)  This  procedure,  typical  for  the  treatment  of  "boundary  layer 
phenomena"  is  always  employed  for  differential  equations  in 
which  the  terms  of  highest  order  sj?e  multiplied  by  small  factors, 
(In  our  case  these  factors  are  yu..  A,  and  3"^  )• 

No  accuracy  would  be  gained  by  trying  to  discuss  the  solutions 
'  X  5  X,  ;  for,  the  conservation  laws  (1) 


for  a  finite  range  x  5  x  -  x,  j 

and  (2)   would  not  be  accurately  valid  iinless  accidentally 

For  the  infinite  range  it  follows  from  regularity  (see  Section  2, 
that  these  derivatives  vanish  at  the  end  points. 
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deduce   that  the  flux  ni  is  positive: 

m  =  mfde  =J    (t  -  iXl  -£)  S(0)dx  >  0  . 

Prom  this  fact  it  follov/s  that  the  reaction  begins  in  the  un- 
burnt  gas  at  x  =  -  oo  and  ends  in  the  burnt  gas  at  x  =  oo  • 

The  values  p  ,  t  ,  9  ,  v  ,  and  p,,  'r.,  0.,  v,  for  unburnt 
ahd  burnt  gas  are  those  that  are  prescribed  at  both  sides  of  the 
discontinuity  front.   These  quantities  are  not  prescribed  arbitra- 
rily, they  are  to  satisfy  the  conservation  laws;  (O)  and 

(1)^  p^  +  mv^  =  p^  +  inv^=  P 

^2)«  Eo*®o-^^^o  =  ^l*®l*i^?  =  ^  ' 

through  which  at  the  same  time  the  values  of  the  constants  P  and 
Q  are  determined*   Here  we  have  set 

E^°^e^)  =  E^  and  E^^^e^^)  =  E^  . 

The  conservation  laws  follow  immediately  for  any  regular 
solution  from  the  differential  equations  (l)'and  (2)',  since  these 
equations  reduce  to  (D^  and  (2)^  for  x  =  oo  and  x  =  -  oo  . 

A  further  assun5)tion  which  we  liispose   on  our  boundary  values 
is 


expressing  that  no  reaction  can  take  place  in  the  unburnt  gas  at 

its  initial  temperature,  and  that  reaction  would  take  place  at 

the  final  terperature  if  unburnt  gas  were  still  left.  ' 

For  the  following  arguments  it  is  convenient  to  eliminate 

v  and  p  by  . 

V  =  mr  and   p  =  r~-^g     , 

1)  Relation  9^  >  9g  would  in5)ly  S  >  0  for  9  =  9^  and  would  hence 
not  be  conpatible  with  equation  (3)'  for  a  regular  solution. 
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and  to  consider    t,    ©   and   £  as   the   only  dependent  variables; 
the  equations  then  become 

(1)  -  yum  g  +  T-^©  +  m\     =  P       , 

^2)  -X^  +  m[E^^^©)   -  |ii?'r^+  Pr]  =  mQ       , 

-  m||  +T"^(l-e)   S(©)    =  0        , 

The   constant  coefficients  m,    P,   Q,    and  the  bo\indary  values 

T    ,     Tn ,   ©   ,   e,  ,    e     =0,    and  e,    =  1  are  subject   to  the  conserva- 
o'        1'      o*      1'      o  '  1  *' 

tion  laws 

(Do  '"o"\-^  ^"^  '^o  =  V\  +  "^^-^1  =  P     » 


v;e  now  consider  the  system  of  differential  equations  (1), 
(2),  (3)  and  ask  whether  or  not  it  possesses  regular  solutions 
assuming  the  prescribed  bovuadary  values. 

2.  Refflilarity  and  determinacy. 

We  shall  introduce  as  "degree  of  regularity"  for  each  end- 
point  a  number  r  such  that  the  manifold  of  solutions  of  the  differ- 
ential equations  (1),  (2),  (3)  which  are  regular  at  that  endpoint 
and  assume  the  prescribed  boundary  values  there  depends  on  r  para- 
meters. 

Let  us  first  consider  the  endpoint  x  =  oo ,  To  determine 
the  degree  of  regularity  r,  at  the  end  x  =  oo  we  shall  introduce 
characteristic  exponents,  oc,  by  the  following  formal  procedure. 
We  expand  the  differential  equations  (1),  (2),  (3)  with  respect 
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to  powers  of  '^ ''\'    ©  -  ©t  »   a^cl  £-£n  =£-1  .   The  terins  of  order 
zero  vanish  by  definition  of  P  and  Q.   The  terms  of  first  order 
constitute  a  system  of  linear  differential  equations  v/ith  constant 
coefficients,  the  "linearized"  differential  equations  for  x  =  oo  . 
Upon  inserting  into  these  equations  multiples  of  an  exponential 
function  e   for  T-Tt*  ©  -  ©-i>  aJ^d  e-l,  we  obtain  three  linear 
equations  for  the  three  coefficients  whose  determinant  is  a  cubic 
equation  for  oc  ,  The  three  roots  of  this  cubic  equation  are  the 
three  characteristic  exponents.  Me   anticipate  the  fact,  shovm 
later  on,  that  for  our  equations  these  characteristic  exponents 
are  always  real.  Suppose  we  have  r^  negative  and  3-r-  positive 
characteristic  exponents,  then  r,  is  the  degree  of  regularity. 
This  fact  is  implied  by  the  well-lmown  theory  of  singular  points 
of  ordinary  differential  equations.  The  behavior  of  the  r, -para- 
metric set  of  regular  solutions  can  be  characterized  by  the  regular 
solutions  of  the  linearized  equations;   the  latter  are  a  linear 
combination  of  exponential  functions  e"''"'^  if  the  characteristic 
roots  are  different;   otherv/lse  terms  like  xe"'^  or  7i  e'^   enter. 
In  case  one  characteristic  exponent  is  zero,  one  cannot  say  off- 
hand what  the  degree  of  regularity  is;  a  special  consideration 
is  needed. 

The  linearized  equations  at  x  =  oo f or  the  quantities 
T-T,,  9-a,  ,  £-1   are  immediately  found  to  be 

.          dC-r-T  ) 
(1)       m/^^  -^^   +  Tj_  ^O^(T-T^)  -  m^(T-tj^)  -r^  10-Q^)=  0, 

(3j.      X,^i^.m|^.AE,(..l).r-X(r-T,)]  =0, 


-  m- 

l^^i-i 

-AE^ 

(£-1) 

..-1, 

m^ 

dx 

.  .- 

1       _ 
-1 

~T© 

i£i  , 

AE  = 

E<°'- 

(3)^  ^^^^T^^    '^r  Si(e-i)      =0, 

where  we  have   set     -5^  =  """^  ^r^^'    ,    AE  =  E^""'-  E^-"-^    ,    and  the 
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and  S,    are  to  be   taken  for  ©  =  0^,    r  =  r^,  £=  1    , 

For  the   characteristic   exponent  oc  one   then  obtains   the 
equation 

1  u.-,  oc  -  m     +  T",    6-, 


-I' 

0 

^- 

m^E^ 

0 

moc 

H-V^Sj 

=  0 


(4)  [mT^06+  S^][mT^  ^/^f-^  '   i^^-^l^  "  ©i\ 


The  first  bracket  has  evidently  one  negative  root.  As  to  the 

2 
second  bracket,  which  we  write  in  the  form  aoc  -  hoc-^   c  we  observe 

that  it  has  real  roots  since  the  discriminant 

2  2. 

b^  -  4ac  =  (m^T^X^-  Q^K^+  m^r^  j^f   +  4  "^  ^^  j'^-^Vl""^^^ 

is  positive. 

Since  a  >  0  we  see  that  the  second  bracket  has  one  positive 
and  one  negative  root  if  c  «:  0,   If  c  =  0,  we  have  b  =  b  = 

(V,-1)Q-,X^  +  -^  ■_  ^  ®iAS  ^  ^  ^^'^  hence  the  bracket  has  one  vanish- 
ing and  one  positive  root.  If  c  >  0  we  have  b  >  b  >  0  and  hence 
the  bracket  has  two  positive  roots. 
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In  case  c  <  0  we  have  r,  =  2  and  in  case  c  >  0  we  have 
r,  =  1.   A  detailed  Investigation  of  the  vector  field  correspond- 
ing to  the  differential  equation  would  show  that  in  case  c  =  0 
a  two  parametric  set  of  regular  solutions  exists;  hence  r^  =  2 
also  in  this  case.   Since  the  condition  c  $  0  is  equivalent  with 
^,Q-  ^  m^r?  we  have  to  distinguish  the  following  two  cases 

Case  (A^)    m^T^  >  V^O^  , 

2„2 


Case  (B^)    ra^^   -\'^i 
The  degree  of  regularity  in  these  cases  is: 


Case  (A,) 


:ase  (B-j^) 


Employing  the  sovmd  speed  c,  =  /  ^t  ®t  ^'^^  ^^^  bxirnt  gas  in 
state  (1)  we  can  write  the  condition  for  cases  (A, )  and  (B, )  in 
the  form 

Case  (A^)    ^i  >  °i   » 

Case  (B^ )    ^l  -  ^1   • 

Thus  the  flow  of  the  burnt  gas  in  state  (1)  is  supersonic 
in  case  (A,)  and  subsonic  or  sonic  in  case  (B,). 

For  the  state  (o)  at  the  end  x  =  -  oowe  obtain  a  similar 
equation  for  o6,  the  only  difference  being  that  S  =  0  in  state 
(o)  since  ©  -  euwas  assumed.  The  equation  for  <x  then  becomes 

(4)-     moc[mt^>^^„<x2-  (m^r^x^-  0^\   +  m^-,^  ^  )=. 
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where  X^,  JJ-^,    V^  refer  to  0  =  9^ ,  'T  =  T^,  &=  0.  The  first 
factor  here  has  the  root  «.=  o.  The  bracket  always  has  one 
positive  root  and  in  addition  one  positive,  vanishing,  or 
negative  root  depending  on  whether  VqOq-  m  ^  >0,=  0,or  <  0, 

;  the  end 
Hence  for 

every  regular  solution  assviming  the  prescribed  boundary  values 
at  X  =  -  CQ  we  have  ©  <  ©gif  x  is  sufficiently  negative,  for 
-  00  <  X  5  X  say.  By  virtue  of  S  =  0  for  ©  <  ©^  equation  (3) 
entails  £=  0  for  -  <X)<  x  -  x  .  The  investigation  of  the  mani- 
fold of  solutions  regular  at  x  =  -  oo  is  thus  reduced  to  the 
investigation  of  the  regular  solutions  of  equation  (1)  and  (2) 
with  E^^^ (©)  =  E^°^(©).  The  characteristic  exponents  for  this 
problem  are  the  roots  of  the  bracket  in  relation  (4)",  Hence 
we  conclude:   If  both  roots  of  the  bracket  are  positive  we  have 

r  =  2,  The  same  is  true  if  one  of  the  roots  is  zero,  as  a 
o 
detailed  investigation  of  the  vector  field  corresponding  to  the 

differential  equation  would  show.   If  one  of  the  roots  is  nega- 
tive, however,  the  other  one  being  positive,  we  have  r^=  1. 
Accordingly,  we  distinguish  the  following  two  cases: 

case  (A^)  m^r^  ^  ^^©^  , 


o 
The  degree  of  regularity  is 


Case  (B^)  m^rj  <  ^^©, 


r  =  2    in  case  (A^)  , 

r  =  1    in  case  (B  )  • 

Employing  the  sound  speed  c  in  the  unb\irnt  gas  we  write 
the  condition  for  cases  (A^)  and  (B^)  in  the  form 

Case  (A^)  v^  >  c^ 

Case  (B^)  v^  <  c^   . 
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Thus,  the  flow  of  the  unburnt  gas  at  x  =  -  oo  Is  supersonic 

or  sonic  there  in  case  A  ,  while  it  is  subsonic  in  case  B  • 
o  o 

We  see  that  four  different  cases  are  to  be  distinguished 
according  to  whether  case  A  or  case  B  obtains  at  x  =  ooor  at 
X  =  -  00  .   It  is  now  very  interesting  that  by  virtue  of  Jouguet ' s 
rule  (see  Supplement  [3],  p.  10)  these  four  cases  just  corres- 
pond to  the  four  cases  of  strong  and  weak  detonations  and  defla- 
grations; the  Chapman- Jouguet  detonation  or  deflagration,  charac- 
terized by  V-,  =  c,  ,  is  here  classed  with  the  strong  detonations 
or  weak  deflagrations  respectively. 

We  now  shall  determine  in  a  formal  way,  simply  by  counting 
the  number  of  parameters,  the  manifold  of  solutions  which  are  re- 
gular at  both  endpoints.   The  set  of  solutions  which  are  regular 
at  X  =  00  is  r, -parametric.  Among  these  solutions  those  are  to  ^ 
be  selected  which  are  regular  at  x  =  -  oo .  Since  all  solutions 
regular  at  x  =  -  cd  form  a  r  -  parametric  set  in  the  three- 
parametric  set  of  all  solutions,  it  is  clear  that  the  condition 
to  be  regular  at  x  =  -  oo  is  expressed  by  3-r^ relations.  Thus 
3-r  conditions  are  imposed  on  the  r,  paj*ameters  characterizing 
the  solutions  regular  at  x  =  oo  ,   One  of  these  parameters  can 
always  be  chosen  arbitrarily  (within  limits),  since  from  every 
solution  satisfying  the  boundary  conditions  one  obtains  a  set 
of  others  by  substituting  x  + const,  for  x.  Thus  3 -r^  conditions 
are  imposed  on  Tt-I  parameters.   If  r-j^-  1>  3-rQ,an  (1"^+  r^-  ^- 
parametric  set  of  solutions  can  be  expected  to  exist.   If 

r,  -1  =  3-r  ,  one  solution  (or  else  a  finite  number  of  them)  can 
1       o  ' 

be  expected  to  exist.   If  r,  -  1  <  3  -  r  more  conditions  are 
imposed  than  parameters  are  available.  These  conditions  will  be' 
satisfied  only  if  the  coefficients  entering  the  differential 
equations  or  the  boundary  values  assume  appropriate  values.  In 
other  words,  4  -  r  -  r,  conditions  are  imposed  on  coefficients 
and  boundary  values.  We  term  the  number 


the  "degree  of  over-determlnacy."  Prom  Jouguet* 3  rule  and  the 
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determination  of  the  values  of  r  and  r,  given  before  we  find 
Strong  detonation,    case  (A^B^),    s  =  0 
Weak  detonation  ,    case  (A^A,),    s  =  1 
Weak  deflagration,    case  (B  B,),    s  =  1 
Strong  deflagration,  case  (B  A^),    s  =  2 

Upon  comparing  this  table  with  the  table  for  the  degree 
of  under-determlnacy  given  In  the  introduction  (p. 6  )  we  realize 
the  fundamental  fact  that  the  degree  of  over-determinacy  re- 
sulting from  the  Internal  mechanism  of  the  reaction  process 
e(:iuals  the  degree  of  under-determlnacy  of  the  external  problem 
flow.   The  number  of  conditions  needed  to  make  the  flov^  problem 
unique  thus  equals  the  number  of  conditions  Imposed  by  the 
mechanism  of  the  reaction  process. 

This  result  may  be  interpreted  as  follows:   All  strong 
detonations  are  possible.  WeaJc  or  Chapman- Jouguet  detonations 
are  only  possible  If  one  of  the  parameters  of  the  process  satis- 
fies one  condition.   As  such  parameter  we  may  consider  the  flux 
m.   As  we  shall  see  later,  weak  detonations  exist  just  for  such 
values  of  the  flux  that  lead  to  Chapman-Jouguet  detonations  ex- 
cept for  sufficiently  high  values  of  the  reaction  rate,  for 
which  a  larger  value  of  the  flux  leads  to  a  possible  weak  deto- 
nation. Weak  deflagrations  are  also  only  possible  if  one  of 
the  parameters,  the  flux  m  say,  satisfies  one  condition.   As 
we  shall  see  later,  wealc  deflagrations  exist  indeed  for  only  a 
particular  value  of  the  fliox.   Strong  deflagrations  should 
exist  only  if  two  conditions  are  satisfied  by  the  parameters. 
As  we  shall  see  later,  strong  deflagrations  do  not  exist  at  all. 

It  must  be  emphasized  that  these  statements  are  so  far 
derived  in  a  purely  formal  manner.  They  are  obtained  by  balanc- 
ing the  number  of  available  parameters  with  the  number  of  condi- 
tions imposed.   A  definite  statement  about  existence  and  unique 
ness  cannot  be  made  on  this  basis.   As  a  matter  of  fact  t-hese 
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arguiients  are  not  sufficient  to  exclude  weak  detonations 
and  strong  deflagrations,   tore  detailed  considerations  are 
needed  for  this  purpose.   In  the  following  Sections  (3  and  4) 
we  shall  first  investigate  certain  limiting  cases  of  weak  de- 
flagrations and  strong  detonations,  and  then  proceed  to  discuss 
in  Section  5  the  problem  of  existence  of  solutions  In  general. 

3«   Weak  deflagrations. 

The  degree  of  over-deterinlnacy  was  found  to  be  s  =  1  for 
v/eak  deflagrations.   As  indicated  above,  such  deflagrations 
can  therefore  be  expected  to  exist  only  if  the  coefficients  of 
the  differential  equations  satisfy  one  condition.   This  condi- 
tion may  be  considered  a  condition  for  the  constant  m,  the  flux, 
or  for  the  burning  velocity  v  =  't'  m.   The  question  arises 
whether  or  not  to  given  values  of  r  ,  ©  ,  and  given  functions 
S(e)  and  E^'^'(e)  there  really  is  a  value  of  the  flux  m  for  which 
the  problem  has  a  solution.   We  shall  deal  with  this  question 
in  Section  5. 

Here  we  shall  consider  only  the  limiting  case  of  a  constant 
pressure  deflagration  in  which  this  question  can  easily  be  an- 
swered.  This  limit  case  results  if  either  heat  conductivity  X 

and  viscosity  u   approach  zero,  if  the  reaction  rate  S  approaches 

1) 
zero   ,  or  if  the  pressure  p^  Increases  indefinitely  while  the 

temperature  9^  remains  fixed.   This  limit  case,  to  which  actual 

situations  frequently  come  very  near,  has  always  been  assumed  for 

numerical  determinations  of  the  flame  speed  (see  [ll],[l2]). 

1)  For  example  by  letting  S   — >  0,  (cf.  footnote  on  p.  12  ). 

The  objection  may  be  raised  that  for  small  values  of  the 
reaction  rate  the  basic  assunption  (pp.  3  ond  10)  cannot  be 
satisfied.   Nevertheless  it  is  necessary  to  consider  the  limit 
S  — >  0  for  the  discussion  of  the  mathematical  implications  of 
the  problem  formulated  in  Section  1. 
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To  describe  the  limiting  "constant  pressiire  problem" 
we  first  introduce  dimensionless  quantities 


P-   ='*S  /p  m^    ,    X  =  S  /p  m^ 
/     oo'  ^o     »       oo'  ^o 


P  =  P/i^Q   ,  Q  =  Q/e^ 


and  Introduce  the  new  variables 


x  =  xS^/mr^ 

The  equations  (1),  (2),  (3)  then  become 

-um^dr/dS  +  ©/'P  +  "^^  =  P 

-XdQ/dx  +  e^^^0^^)/9q  -  ifi^'f^  +  PT=  Q 
-  de/dx  +r-  1(1  -e)  S(0^e)/©Q       =  0  , 

This  system  of  differential  equations  depends  on  the 
parameters  m,  X,  ^u,  P,  Q,,  and  ©  .  We  obtain  the  limiting  pro- 
blem by  considering  a  set  of  such  systems  for  v;hich  m  approaches 
zero  while  the  other  parameters  are  kept  fixed.   In  other  words: 
The  differential  equations  of  the  limiting  problem  are  simply  ob- 
tained  by  omitting  the  three  terms  involving  the  factor  m  . 

The  fact  that  X  and  "Jx   are  fixed  while  m  approaches  zero 
evidently  implies  ^S^/p^— >0  and  /xS^/p^ -*»  0  ,  The  limiting 
equations  will  therefore  represent  a  good  approximation  if 
As  /p^jyu-S  /p  ,  and  m  are  small.   In  that  case  we  can  re- 
introduce the  original  quantities.   Thus  we  obtain  the  equations 
of  the  constant  pressure  problem  in  the  following  form: 
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(T) 


(2)  -y^  +  m  [E^^Uq)   +0]  =  m  Q 


d 


(5)  -31  4=  +  PQ'-"-  (1  -  ^  )S(9)  =  0 


where  we  have  eliminated  t"  f roin  (2)  and  (3),  using  Cl). 

Relation  (l)  expresses  the  fact  that  under  the  condi- 
tions of  the  limiting  problem  the  pressure  does  not  vary  across 
the  flame  front:  we  have  constant  pressure  combustion.   Fur- 
ther we  see  that  the  term  representing  the  kinetic  energy  has 
dropped  out  from  (2). 

Tlie  boundary  conditions  are  9  =  ©,,  £  =  1  at  x  =  od  , 
and  £=Oatx=-oo.   The  constant  Q  Is  given  by  Q  =  E.  +9,  . 
The  relation  ^^  +  9^  =  Q,  which  holds  for  regular  solutions  deter- 
mines the  value  9  =  9  atx  =  -oo. 
o 

To  investigate  possible  solutions  of  this  problem  one  may 
consider  9  as  Independent  variable  running  from  9  to  9,  and 
combine  the  equations  to 

(^)  dL      2A  S(9)  1  -  £ 

^'^  ^   ^-^E,  .9,  -E^^^9)  -9 

It  is  immediately  seen  that  this  equation  has  a  saddle-sliagularity 

at  the  point  9  =  9^,  £  =  1.   There  is  just  one  solution  curve 

that  enters  this  point  from  the  region  £<  1,  9  <  9,.   If  this  curve 

Is  followed  backwards  it  will  enter  the  axis  £  =  0  at  a  point 

with  9  >  ©3  provided  p  A/m^  is  sufficiently  large;  If  p  A/m^  is 

sufficiently  small  the  curve  will  enter  the  line  9=9  at  a 

point  with  E  >  0,      It  is  thus  clear  that  there  is  just  one  value 

of  p  A/m^  for  which  the  curve  enters  the  line  9=9  with  £  =  0. 

Since  d£/d9  =  0  for  9  <  9^ ,  the  curve  will  enter  the  line  9=9 
'  —  3  o 

also  with  the  value  c  =  0.   That  for  this  solution  the  quantity 
X   approaches  oo  as  9  -*>  9,  and  -  cd  as  9  — >9  is  Immediately  seen 
from  (2)  and  (3). 
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Thus  it  is  shown  that  in  the  limit,  as  the  quantities 
Xe  /^  anduS  /p  approach  zero,  a  deflagration  process  exists 
with  a  welldef ined  flux  m  and  flame  speed  v  .  As  will  be 
shown  later  (at  the  end  of  Section  5),  the  same  is  true  for 
small  values  of  these  quantities,  in  other  words,  for  sufficiently 
small  V  alues  of  Xs  Vp  and  /^Sjn/Po  ^  deflagration  exists  which 
is  nearly  a  constant  pressure  combustion  and  which  is  character- 
ized bj  a  uniquely  defined  burning;  speed  v^ 

Computations  of  the  burning  speed  can  be  carried  out  by 
solving  equations  (2)  and  (3)  or  (5)  through  interactions,  (see 
[10],  and  [11],  in  the  latter  report  the  term  d(  S  dg/dx)dx  was 
added  to  equation  (3)  in  order  to  take  diffusion  Into  account); 
it  was  found  in  the  quoted  reports  that  for  the  actual  situations 
considered  the  approximation  involved  in  (2)  and  (3)  was  rather 
accurate  since  the  omitted  terms  turned  out  to  be  very  small 
for  the  calculated  solution. 

It  is  interesting  that  essentially  only  the  combination 
pA/m  or  the  diraensionless  combination 

(6)  A/P^^OD  =  ^^>o^/^>^Sqd^  ^ 

enters  the  constant  pressure  problem,  as  seen  from  equation  (5). 
(More  precisely,  the  problem  depends  only  on  the  dimenslonless 
quantity  (6)  in  addition  to  the  functions  S(9)/S^  and  E^^'(©)/©q; 
note  E  +  ®Q=  E,  +  ft,  . )   The  quantity  (6)  approaches  a  finite 
limit  value  as  AS  /p  approaches  zero.   Therefore,  for  small 


values  of  As  /p  ,  the  flux  m  is  proportional  to  /p^X^/^q  and 

the  burning  speed  v^  ^o  A^oo®o/Po'   •'"^  ^^  particular  X  ,^  ,  sxid 

the  reduced  initial  temperature  0^  are  kept  fixed,  the  flux 

increases  like  /s/p  and  the  burning  speed  decreases  like 

j/sZTp"  as  p  increases.   If  S^  v/ero  independent  of  p  or  increased 

with  p  of  less  than  first  order,  the  latter  result  would  mean 

o 
that  the  burning:  speed  sgould  decrease  with  increasing  pressure. 

This  does  not  seem  to  be  confirmed  experiment ally D .  It  is  thus 


1)  See  e.g.  [9]  p.  146,  Table  27< 
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indicated  that  S^  increases 
oo 


ing  expression  (6),  since  the  approximation  made  is  no  longer 
valid  if  PqA^  Is  small.  Whether  or  not  deflagrations  are 
then  possible  will  be  discussed  later  on,  in  Section  5, 


4,  Strong  detonations  and  Chapman- Joup;uet  detonations. 

We  next  consider  strong  detonations,  including  Chapman- 
Jouguet  detonation;  according  to  Jouguet's  rule  they  have  it 
in  common  with  weak  deflagrations  that  the  flow  of  the  burnt 
gas  is  subsonic.  The  flow  of  the  unb\irnt  gas, , however,  is 
supersonic  (or  sonic).  The  degree  of  over-determlnacy  was  found 
to  be  3  =0,  Hence  no  conditions  need  be  imposed  on  the  data 
other  than  the  inequality  m  r   <  ^i®i»  °^   ^i  -  '^i*  V/hether  or 
not  it  is  true  under  this  condition  that  a  unique  solution 
exists  will  be  discussed  in  the  next  Section,  5, 

In  the  present  Section  we  shall  consider  a  limiting  case 
obtained  by  letting  yb»-s/p  and  X  sVp  approach  zero.  We  imagine, 
in  particvilar  that  viscosity  yu  and  heat  conductivity  X  approach 
zero  while  S  and  p  are  kept  fixed.  This  limiting  problem  pre- 
sents an  analogy  to  the  limiting  problem  considered  in  Section 
3;  the  limiting  behavior  of  the  present  limiting  problem  is, 
however,  quite  different  from  that  encountered  in  Section  3, 

Consider  a  sequence  of  solutions.  There  are  two  possi- 
bilities:  either  the  terras /a-?-  andX-?—  in  equations  (2)  drop 
out  in  the  limit,  or  -?-  a^^dg*  become  infinite.   For  those  values 
of  X  for  which  the  first  case  occurs,  the  equations 

(1)^  t"^©  +  m^T    =  P  , 

(2)  E^^^e)-  imV  +P'r  =  Q  , 

00  ii 

(3)  -  °^li  +  '^"^^1  -')  ^^®^   =  °  ' 
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are  satisfied  in  the  limit.   If  the  second  case  occtirs  at  a 
place  X,  a  discontinuity  of  t  and  0  occiirs  in  the  limit  while 
£,  remains  continuous.  Such  a  discontinuity  would  simply  be 
a  shock  not  involving  a  reaction, 

A  more  detailed  investigation  (see  Section  5)  of  the 
limit  process  ^S^^p  ,  /"-S  /p  -*'  0  will  yield  that  a  strong 
detonation  can  in  the  limit  be  described  as  a  shock,  not  involv- 
ing a  reaction,  immediately  followed  by  a  reaction  process  gov- 
erned by  equations  (!)__>  ^^^m*  ^^^*  ^^^  confirms  the  accepted 
idea  about  a  detonation  process,  which  was  formiilated  more  speci- 
fically by  G.  I.  Taylor  and  v.  Neximaim  (see  [1],  e.g.). 

The  question  arises  whether  to  given  values  of  T  ,  ©  ,  r^, 
-1  —1  o   o   i 

ft.  ,  and  m  with  T   c  <  m  5  t^  c,  satisfying  the  conservation 

laws,  a  transition  process  exists  which  consists  of  a  shock 

followed  by  a  reaction.  We  assume  that  the  temperatvire  to  vrtiich 

the  shock  raises  the  unb\irnt  gas  is  above  safety  temperatvire. 

Otherwise  the  reaction  process  would  simply  be  a  deflagration 

process,  which  cannot  exist  in  the  present  limiting  case   . 

Denoting  the  quantities  past  the  shock  front  by  an  asterisk  we 

require 

%    <%    <   »^.    • 

To  find  out  whether  the  equations  (1)  ,  (2)  ,  (3)  possess  a 
solution  we  determine  T  and  ©  as  functions  of  e  from  (l)-^  and 
(2)  ,  and  insert  in  (3),  This  is  possible  if  the  Jacobian 

(m*^  -  T  '^e)  -  T-'^Q 


d© 
of  (1)   and  (2)   does  not  vanish.   Introducing  quantities  V^ 


^(e) 


dE^^V<a©=  i/(/^^-  1) 


1)   Expressing   r  and  ©  through  e  by    (1)„  and   (2)^    equation   (3) 
becomes    a  differential  equation  f or  £  ,^hlch  yi^ds   de/dx  =  0 
for   £  =  0  since  S(eb)   =  0  for  ©b    5  ©s .     The  sole   solution  of 
this    differential   equation  vanishing  for  x  =  -  oois   therefore 
d    =  0    . 
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we  find  p     p  „ 

J  =  (m'^  -  T-2c^)/(^-  1)    , 

where  T ,  c,  and  V  depend  on  e  .  Since  tne  flow  is  subsonic 
in  the  state  (*)  past  the  shock  we  have  J^  <  0.   In  the  neigh- 
borhood of  the  state  (*)  we  may  therefore  express  r  and  9     in 
terms  of  £  and  the  solution  of  the  equation  resulting  from  (3) 
is  uniquely  determined  through  the  initial  condition  e  =  0  at 
X  =  0,  say,  where  we  may  place  the  shock  front.  The  question 
then  is  whether  on  contin\iing  the  solution  one  would  obtain  a 
value  for  which  J  changes  sign.  The  final  state  is  also 
subsonic,  J,  <  0,  and  the  transitions  from  the  state  (■«•)  to  any 
of  the  intermediate  states  between  (•j:-)  and  (1)  correspond,  as 
regards  the  conservation  laws,  to  a  set  of  weak  deflagrations. 
All  the  intermediate  states  are  thus  subsonic.  Consequently, 
J  remains  negative  throughout   , 

The  differential  equation  for  e   resvilting  from  (3) 
possesses,  therefore,  a  solution  with  e  =  0  for  x  =  0,  Since 
for  this  solution,  de/dx  -*•  0  as  £  -*>  1,  the  final  state  (1) 
is  approached  as  x  -»►  oo  .  Thus  it  is  seen  that  a  limiting  type 
detonation  consisting  of  a  shock  followed  by  a  reaction  is  poss- 
ible for  arbitrary  values  of  the  flux  m,  satisfying  'v'j^c     < 

The  reaction  process  in  a  strong  (limit  type)  detonation 
following  the  shock  has  it  in  common  with  a  weak  deflagration 

1)  Whether  or  not  dE^^'/de   changes  sign  during  this  process  is 
Immaterial,  cf,  however,  v,  Neumann*s  report  [1],  in  which  transi- 
tion processes  are  discussed  on  the  basis  of  equations  (1)   and  (2), 

2)  For  numerical  determination  of  such  processes  and  various  de- 
tailed discussions  see  the  reports  by  Eyring  and  his  collaborators, 
[6] (where  an  excessively  high  reaction  rate  was  assumed)  and  [7], 
(The  case  shovm  in  Pig.  10a,  p,  44  in  [7]  and  labeled  "steady 
deflagration"  there  is  a  "weak  detonation"  in  our  t erminology.) 
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that  unbiirnt  gas  in  a  subsonic  state  is  transformed  into 
■bxirnt  gas  in  a  subsonic  state.   The  two  processes  differ, 
however,  in  other  respects.   The  initial  temperature  in  a 
deflagration  is  belov/  safety  temperature  while  the  reaction 
in  a  detonation  begins  v;ith  a  higher  temperature.  Also  in 
a  deflagration  the  rates  of  change  dV<ix»  dQ/dx,  de/dx  are 
zero  initially  while  there  is  no  such  restriction  for  the 
reaction  process  following  a  shock;   for  small  values  of  X 
and  yu-  the  rates  of  change  tmdergo  such  great  changes  in  the 
narrow  shock  zone  that  adjustment  to  any  values  of  these  rates 
at  the  beginning  of  the  reaction  zone  is  possible.   It  can 
thus  be  \inderstood  that  for  the  reaction  process  in  a  detona- 
tion no  definite  value  of  the  fliix  is  required  as  for  proper 
deflagrations. 

5,  Discussions  of  possible  processes  with  the  aid  of  the  vector 
field. 

It  is  very  helpful  for  the  further  discussion  to  refer 
to  the  vector  field  in  the  (T,  Q,   e.) -space  generated  by  the 
differential  equations  (1),  (2),  and  (3),   It  is  convenient  for 
our  discussion  to  consider  as  positive  the  direction  of  decreas- 
ing x;   when  v/e  speak  of  following  a  solution  curve  sf  we  always 
imply  that  we  travel  in  direction  of  decreasing  x.  All  vectors 
of  the  field  point  to  smaller  values  of  e.  except  on  the  surfaces 
6=1  and  in  the  region  9^0,  v;here  they  lie  in  the  svirf  aces 
£  =  const.  Every  solution  curve  that  hits  the  svirface  9  =  6g 
will  stay  on  the  plane  e  =  const,  from  there  on, 

Ptirther  important  properties  of  this  vector  field  are 

2  2 
these:   dT=  0  on  the  cylinder  ©  +  m  T  =  Prand  d©  =  0  on  the 


surface  ^  given  by  the  equation 

E^^he)  =  |mV  -  PT  +  ^  , 

Points  of  Intersections  of  the  siorface  ^  with  the  cylin- 
2  2 
der  ©  +  m  r^  =  PT  and  the  plane  e.  =  const,  will  be  denoted  by  A  g. 
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and  Bg   •   The  conservation  lavrs    (D   and  (2)   are  obviously 
satisfied  at  such  points.   Any  two  such  points  on  the  plane 
fi  =  const,  with  9  >  0  andT>  0  evidently  represent  the  two 
states  on  both  sides  of  a  possible  shock  transition  in  the 
gas  mixture  characterized  by  the  value  of  e   considered.   There- 
fore, one  of  these  two  points,  Ag  »    corresponds  to  a  super- 
sonic, the  other,  Bg  ,  to  a  subsonic  flovif.  We  assume  the  values 
of  the  constants  P  and  Q  such  that  the  surface  $    intersects 
the  initial  plane  £  =  0  at  two  points  A  and  B  with  ©  >  0,'r>  0, 
Prom  the  discussion  on  p.  23  it  follov;s  that  for  every  £  >  0  two 
points  of  intersection  Ag  and  Be  exist  as  long  as  Ae  or  B^ 
do  not  become  sonic,  or,  what  is  equivalent,  do  not  coalesce. 
We  assume  that  this  is  not  the  case  for  0  <£<  1;  this  assumption 
is  primarily  a  condition  on  the  flux  m.  We  also  assume  that 
0  >  0,  T  >  0  at  these  points  for  0  «  £  5  1.  We  then  have  two 
curves  >A.   and  ©   of  points  Ag  and  Bg  along  which  0  and  T  are 
continuous  fiinctions  of  the  paj-ameter  £ ,  (The  last  assximption 
made  here  is  somev/hat  stronger  than  necessary.  For  the  discussion 
of  weak  deflagrations,  for  example,  only  the  existence  of  the 
curve  0  is  needed.   Incidentally,  the  existence  of  the  super- 
sonic state  Ag  v/ith  0  >  0,-^>  0  alv/ays  implies  the  existence 
of  the  subsonic  state  Bg  >   hut  the  converse  is  not  true.) 

At  the  points  A^,  Be  we  have  dT  =  d©  =  0,  de/dx  >  0  ex- 
cept for  £  =  1  or  ©  5  ©„;  hence  the  field  vector  points  in  the 
negative  £-  direction  at  these  points.  The  projections  of  the 
field  vectors  on  the  planes  £  =  const,  have  singularities  at 
the  points  Ag  and  Bg   ,  At  a  point  Ag.  the  projected  field  has 
a  nodal  point  and  the  solution  CTirves  of  the  projected  field 
lead  from  the  neighborhood  of  this  point  into  it.  At  the  point 
Bg  the  projected  field  has  a  saddle-singularity;  (see  Figures 
2  to  5),  At  the  points  Ag  and  B^  on  the  siirface  £  =  1  and 
in  the  region  ©  5  ©^  the  three-dimensional  vector  field  has 
singularities. 
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(Integral   cvirves   oT 
on  the   plane   € 


Fig.  2 

Detonation 

the  projected  vector  field 
1  for  a  detonation) 
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Pig.  3 
Detonation 

(Integral  ciorves  of  the  projected  vector  field  on 
the  plane  E   =  const,  for  a  detonation) 
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Pis.  4 

Deflagration 

(Integral  curves  of  the  projected  vector  field  on 
the  plane  C  =  1  for  a  deflagration) 
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Pig.  5 
Deflagration 

(Integral  curves  of  the  projected  vector  field  on 
the  plane  £    =  const,  for  a  deflagration) 
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If  we  are  interested  in  a  deflagration  we  must  assume 
that  0  <  0^  for  the  initial  point  B  ,  (see  p.27);   if  we 
are  interested  in  a  detonation  we  must  assume  , (seep, 28  and  Pl^s  .2,  ^ 


Y/e  require  somev/hat  more  for  detonations,  viz,  that  ©  >  8^  on 
the  whole  line  ©  .  Vi/e  first  investigate  which  strong  deflaf;ra- 
tions  and  weal^  detonations  are  possible.   Both  processes  have 
it  in  common  that  the  flow  in  the  biirnt  gas  is  supersonic. 
The  state  (1)  thus  belongs  to  the  case  (A),  and  corresponds  to 
a  point  A-,  .   As  was  shovm  earlier  there  exists  in  case  (A-,)  only 
a  one-parametric  set  of  solutions  v/hich  are  regular  at  x  =  od  , 
and  the  values  of  the  parameter  may  be  chosen  arbitrarily  pro- 
vided it  is  so  chosen  that  e  decreases  as  x  decreases.  Thus 
there  exists  only  one  solution  cui've,  ^  ,  starting  at  A-,  which 
could  represent  one  of  the  processes  mentioned.   If  this  solu- 
tion curve  {f  reaches  the  point  B  it  represents  a  strong  defla- 
gration, if  it  reaches  the  point  A  it  represents  a  wealc  detona- 
tion. 

In  the  follov/lng  we  shall  consider  problems  differing  in 
the  reaction  rate  S,  all  other  parameters  being  unaltered;  we 
may,  for  example,  assume  the  factor  S   (cf.l  on  p.  12)  to  vary 
from  0  to  00  .   The  ci.u'ves  J\r   and  0  are  evidently  independent 
of  S;   the  curve  ^  ,   hov/ever,  depends  on  S.   If  we  want  to 
emphasize  this  dependence  we  write   (S), 

If  the  reaction  rate  S  is  very  high  compared  with  r/n 
and  Pq/u-  ,    all  vectors  point  nearly  in  negative  &  -direction 
except  near  6=1  and  9  ^  ©^ ,   In  this  case,  therefore  any 
solution  c\arve  is  approximately  a  straight  line  in  the  negative 
e  -direction  until  it  meets  the  surface  ©  =  ©  ,  Therefore, 
the  curve  <^  that  begins  at  the  point  A,,  i.  e.  at  c  =  1,  ©  =  Q-j^  j 
T  =  T  ,  ends  up  on  6=0  nearly  with  the  values  ©  =  ©,  ,  T  =  t 
hence  with  nearly  9  =  ©-,  >  9  ,   Thus  we  see  that  f^    ends  up  on 
6  =  0  at  a  point  with  9  >  9  if  S  is  sufficiently  high.   At 
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such  a  point  dg/dx  >  0  and  hence  there  Is  no  continuation 
of  the  curve  on  the  plane  e  =  0  and  hence  none  of  the  desired 
initial  states  is  reached. 

Let  us  consider  the  opposite  extreme  that  the  reaction 
rate  S  is  very  small;  then  the  direction  of  the  field  vector 
lies  everyv/here  nearly  in  the  plane  &  =  const,  except  near  the 
curves  ^ and  iu  •  Consider  any  "cylindrical"  neighborhood  of 
the  line  Jt    and  exclude  from  it  an  arbitrarily  small  neighbor- 
hood of  the  point  A,  ,   If  S  is  small  enour;h  then  clearly  the 
vector  field  on  the  lateral  surface  of  the  "cylinder"  points 
into  its  interior.   The  cxirve  ^,  beginning  at  A,,  can  there- 
fore never  lead  far  away  from  the  curve o/r;   for,  as  soon  as 
moved  away  from  (^  ,    the  rate  of  change  de/dx  would  become  much 
smaller  than  |dT/dx|  -f-    |d9/dx|  and  therefore  the  curve  ^   would 
again  be  dravm  nearer  into  the  neighborhood  of  v>^.  Consequently , 
the  curve  d^    meets  the  surface  ©  =  ©  not  far  from  the  inter- 
section A  of  this  surface  with  the  c\arve  Jj ,   remains  from  there 
on  the  plane  £  =  const  and  soon  enters  the  point  A  £  , 

In  general  the  course  of  the  curve  ^  can  be  delimited  as 
follows:   Let  (%,©£)  be  the  coordinates  of  the  point  Ag,  ,  From 
the  assumption  9E^^/3e>  0  made  earlier  (p. 20)  it  follows  that  T^ 
increases  and  Qc  decreases  v;ith  decreasing  £.,1)  The  statement 


1)  Prom  (1)  and  (2)  v/e  have 


y,- 1     9 


3f^^^         -1 


'£   «»£  "^  "e    » 


d©£  =  -  (m^-Tg  - 'T'r'''Q£)dTe 

-1     2 

where  P  =  %  ©^  +  m  %  has  been  used.  Eliminating  d©£  one  has 

(  ^^  -  1)^^-  d£  =  (m^r^  -  %r'\^   )dr^ 
whence  the  statement  follorjs  since  the  state  k^     is  supersonic 
and  hence  m  T^  "  ^f  ®£  • 


then  is   that   the   cvirve    (^   remains   in  the   cell 

T^  <  T  <  r^,   ©^   <  Q  <  0^,   0  <  e:  <  1     . 

This  follows  from  the  fact  that  the  field  vector  at  the  boundary 
of  this  cell  points  into  its  interior  except  on  <£.  =  0, 

This  statement  implies  that  the  curve  $  never  ends  up 
at  the  point  B  ,  since  B  is  evidently  not  contained  in  the  cell 
just  described.   In  other  words,  strong  deflap:rations  are  imposs- 
ible. 

Let  us  denote  the  points  at  which  the  curve  ^  enters 
either  the  surface  £  =  0  or  the  surface  ©  =  ©«  as  "terminal  points" 
T  w  .  They  form  a  "terminal  line"  tf  which  connects  the  point  'r  =  'r^ 
0  =  0^,  f  =  0  (for  infinite  S)  with  the  point  A^(for  S  =  0) . 

The  terminal  point  T^  depends  continuously  '   on  the  reac- 
tion rate  S,  Therefore,  if  we  let  S  vary  from  cd  to  0,  the  terminal 

point  will  move  continuously  along  the  terminal  line  from  the 
plane  £  =  0  into  the  surface  0  =  ©  »  Consequently  there  is  a 
particular  value  S  of  S  for  which  the  terminal  point  lies  on  the 
intersection  of  the  plane  £  =  0  and  the  surface  9  =  9^ .  We  then 
speak  of  the  "extreme"  situation  and  the  "extreme"  terminal  point 
Tt  Otherwise,  if  the  terminal  point  lies  on  the  surface  9  =  0^, 
we  speak  of  a  "normal"  situation,  if  it  lies  on  £.  =  0,  of  an 
"abnormal"  one. 

In  the  extreme  situation  a  weak  detonation  exists;  for, 

since  9  =  9  at  the  extreme  terminal  point,  the  solution  curve 

starting  at  this  point  remains  on  the  surface  e  =  0.  That  this 

solution  curve  ends  up  at  A  follows  from  the  fact  that  the  field 
•^     o 

has  an  attractive  singularity  at  A^  and  that  the  terminal  point 
T*  =  (T",  9*)  lies  in  the  cell  T^  ^  ^'-^  \,   p^  S  p»  :C  p^  and  hence 

It  is  clear 

1)  This  follows  immediately  from  d9/dx  7^0  on  9  =9,£>0  and 
dd/dx  ?^  0  on  e  =  0,  9  >  9  ;  the  continuity  at  9  =  ^  ,  £  =  0  can 
be  shovm  by  a  more  refined  consideration. 
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that  a  solution  cxirve  can  reach  A  only  if  its  terminal  point 
lies  on  £  =  0  and  ©  =  ©  .  Hence,  a  weak  detonation  exists  only 
in  the  extreme  situation,  i.e.  if  the  reaction  rate  S  assumes 
a  particular  high  value  S  , 

The  discussion  of  the  terminal  line  c7  will  also  be  use- 
ful for  the  investigation  of  strong  detonations.  Before  we  enter 
the  discussion  of  strong  detonations  and  weak  deflagrations  we 
must  investigate  the  possible  continuations  of  solution  curves 
after  they  have  entered  the  region  ©50. 

In  the  region  ©  5  0  the  actual  vector  field  agrees  with 
the  projected  vector  fiold.  Suppose  the  point  Bg  =  ('^♦©e  ,  &) 
on  the  curve  ©  lies  below  ©  =  ©  ,  i.e.  suppose  ©^  <  ©3*  Then 
the  vector  field  has  a  saddle  singularity  at  Bg  ,  (see  p,3l)« 
Consequently,  two  solution  curves,  oOt  and  eO/,  leave  B^  (see  Pig.S); 
and  two  solution  curves,  ^^  and  y-/" y   enter  B^ ;  (the  curves  with 
©  >  ©£  are  aO/   and  oD^-  ,  the  one  with  ©  <  ©  g  are  D^     and  ^J  ) . 
We  see  from  the  vector  field  (Fig.  3  and  Pig.  5)  that  £}^    ends 
up  at  the  point  Af  on  J\  , 

A  solution  curve  tP   can  enter  the  region  ©  ^  ©  only  at 
a  point  where  d©/dx  5  0,  hence  only  on  the  section  ^  of  the  sur- 
face ©  =  ©  cut  out  by  the  surface  <^(see  p.  30  ).  The  intersection 
of  X  with  a  plane  £  =  const,  will  be  denoted  by  3C^,    (see  Fig,  3 
and  5).  We  consider  the  continuation  of  the  solution  curve  -c/ 
after  its  entry  into  ©  <  ©  on  the  segment  X^,     There  are  two 
cases.  Either,  the  point  B^  =  ('^»©£  >  ^)  lies  above  the  plane 
©  =  ©  ,  i.e.  ©g>  ©  ,  Then  the  continuation  of  -c^  from  X^  on, 
ends  up  in  the  point  A^,  as  seen  from  the  vector  field  (Fig,  3). 
Or,  the  point  B^  lies  in  the  region  ©  <  ©  (see  Pig.  5),  Then 
the  curve  y+  entering  B^  with  decreasing  ©  intersects  the  seg- 
ment ^  in  a  point  Gg  •   If  t7  enters  <^  on  one  side  of  G£(on  the 
side  with  largerr),  its  continuation  ends  up  at  A^  as  before. 
If  Centers  o^^  on  the  other  side  its  continuation  ends  up  on 
T  =  0  Tonless  it  has  left  the  region  ©  5  ©  before  reaching  r=  0, 
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If  -e/  enters  (Z^^  at  G^  the  continuation  of  -^  leads  along  i;^  to 
B£   and  from  there  on  -^  can  be  either  continued  along J0^  up 
to  A£  or  alongoO^  up  to  t  =  0  unless  J3+  leaves  ©  -  ©^  before 
reaching  t'  =  0. 

We  now  enter  the  discussion  of  strong;  detonations.  Strong 
and  Chapman- Jouguet  detonations  imply  subsonic  or  sonic  flow  in 
the  bvirnt  gas  at  x  =  oo  .  Hence  case  (B, )  obtains,  or  the  state 
(1)  corresponds  to  a  point  B^^.  Consequently,  according  to  the 
statements  made  earlier,  (see  Section  3,  p.  19)  there  is  a  two- 
parametric  set  of  solutions  regular  at  x  =  oo  .  Hence  there  is 
a  one-parametric  set  of  solution  ctirves  £^  leaving  the  point  B^. 
We  want  to  show  that  among  the  curves  £f   there  Is  a  one  n\   vrtiich 
ends  up  at  A  ,  and  thus  represents  a  strong  (or  Chapman- Jouguet) 
detonation. 

The  set  of  curves  <D   is  limited  by  two  curves  remaining 
on  the  plane  6=1,   One  of  these  two  curves,  ,0^  ,  leads  to  - 
larger  values  of  ©,  and  ends  up  on  the  plane  T  =  0  at  a  point  Ej^ 
(see  Fig.  2).  The  other  curve, jQj^,  leads  to  smaller  values  of 
9  and  ends  up  in  the  point  A, ;  It  represents  the  possible  shock 
transition  from  state  A^  to  state  Bj_ .  The  set  of  curves  J^   may 
be  characterized  by  a  parameter  /3  ,  which  for  /3  =  0  yields  U^ 
and  for  /3  =  1  yields  jQi^". 

On  the  basis  of  the  assumption  made  earlier  (p. 32)  that 
the  cTirve  Q>  lies  above  the  plane  e  =  Q^  we  see  from  the  remarks 
made  before  that  all  curves  D   that  enter  the  plane  ©  =  ©g  end 
up  on  yhi   none  therefore  meet  S^   below  0  =  ©^  or  end  up  on  ©  =  0, 
Therefore  the  curves  ^   end  up  either  on  'T'=  0,  <£-  =  0,  or  on  JV, 
The  so  defined  end  points  will  be  calleci  "ultimate"  points  and 
denoted  by  U^  (S)  or  U^  (S), 

We  Investigate  whether  or  not  the  point  Ug  depends  con-^ 
tinuously  on  the  parameter  ^   .   The  continuity  of  U^  could  be 
interrupted  only  in  three  cases.  The  first  case  would  be  that 
U/3were  the  ultimate  point  of  a  curve  J0  which  passes  through  a 
saddle  shgularity  of  the  differential  equation.   This  possibility 
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Is  excluded;  for  there  is  no  saddle  singularity  except  on 
6=1   since  fi> was  assumed  to  lie  above  0  =  ©  ,  The  second 
case  wovild  be  that  dr/dx  =  0  at  a  point  U|3  on  Ts  0;   this 
Is  excluded  for  0  >  0  as  seen  from  the  differential  equation 
(1).  The  third  case  would  be  that  dd/dx  =  0  at  a  point  U^  on 
e  =  0,  This  case  can  arise  at  points  on  £  =  0  with  ©  =  0  . 

Let  |3  ^  be  a  value  of  /3    such  that  the  coordinates  & 
and  ©  of  TJa     approach  0  and  ©  respectively  as  /3  t|3^j.  (i.e.  /3 
approaches  /3„  from  below),   A  more  detailed  investigation  would 
show  that  then^'''^  =  JD^  ^^   ends  up  at  a  point  U"'*  =  U^^  with  &=•  0 
and  ©  =  ©g. 

If,  for  values  of  /3  somewhat  greater  than  /3^,   the  c\irve 
cD/3  ends  up  on  £  =  0  with  ©  >  e^,  the  point  U^s  is  continuous 
at  y3^.  If,  however,  such  a  curve  jD^does   not  reach  £  =  0  it 
will  meet  0  =  ©  for  a  value  .f  >  0  and  end  up  on  It.  Thus 
Up  =  A   for  such  values  of  p,  and  as  pip*  (i.e.  as  p  approaches 
f'    from  above),  U„  -*■  A  .■  Thus  Ug  is  not  continuous  at  p  =  p.!;.. 

These  considerations  make  It  evident  that  the  ultimate 
points  U/3  vary  continuously  from  the  position  U^  =  H-j^  on,  until 
either  /3  reaches  a  value  /S^^^   for  which  U"  is  on  £  =  0,  ©  =  ©g 
or  until  for  /3  =  1  a  point  U,  on  £  =  0  with  ©  >  ©^  is  reached 
without  passing  through  such  a  point  U„_.  We  maintain  that  the 
latter  possibility  cannot  arise  in  the  normal  situation  S  <  S  . 

To  this  end  we  investigate  the  curves  J5  for  values  of 
/3  slightly  less  than  1.  They  remain  near  the  curve  <i?^  until 
they  come  near  to  the  point  A, ;  from  then  on  they  will  remain 
near  the  ctirve  <^  ,  which  leaves  the  point  A^.   In  the  limit, 
/3  -c>  1,  we  have  a  cvirve  jO  which  coincides  withc^^"  up  to  the 
point  A,  and  then  coincides  with  the  cxirve^.  The  ultimate  point 
U,,  therefore,  lies  on  e  =  0  with  ©  <  ©g  only  if  this  is  the 
case  for  the  terminal  point  of  ^ ,  viz.  in  the  "abnormal"  situa- 
tion.  In  the  normal  situation  no  terminal  point  on  (S  =  0  with 
©  >  ©  exists.  Consequently,  in  the  normal  situation  there  exists 
an  \iltimate  point  TJ^  on  &=  0  with  ©  =  ©  , 
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The  curve  o(?  =  <J0„  whose  ultimate  point  is  U  remains 
on  the  plane  e  =  0  after  having  passed  through  U^,  and  enters 
the  point  Aq.  It  is  evident  that  the  curve J&^  represents  a 
strong  detonation.  Thus  we  have  established  that  in  normal 
situations  a  strong  detonation  is  possible.  It  is  not  so  easily 
seen  whether  or  not  strong  detonations  are  possible  in  the  ab- 
normal case  since  it  is  not  obvious  whether  or  not  there  are 
ultimate  points  on  e  =  0  with  0=0  from  which  the  continuation 

leads  into  A^. 
o 

Suppose  the  value  of  the  reaction  rate  S  is  such  that  we 
are  in  a  normal  situation  but  near  to  the  extreme  situation.  Then 
the  ultimate  point  U*  is  near  to  the  terminal  point  Tg  of  the 
curve  ip.     Hence  the  CTirve  «5  ending  up  at  U"''"  will  first  lie  near 
to  the  curve  .or  leading  from  B,  to  A,  and  then  near  to  the  curve 
^   leading  from  A,  to  T»  •  Thus  we  see:   in  the  extreme  case, 
in  Tdiilch  the  point  T^  falls  on  T^  ,  the  strong  detonation  is 
represented  by  a  curve  consisting  of  the  c\irveJ?r  in  the  plane 
6=0  from  B,  to  A,,  then  by  (p    leading  from  A,  to  T  ^  and  finally 
by  a  section  leading  from  Tg  to  A  ,  In  other  words:  in  the  ex- 
treme situation  not  only  a  weak  detonation  exists  but  also  a 
strong  detonation  which  consists  of  the  weak  detonation  followed 
by  a  shock.  In  a  normal  but  nearly  extreme  situation  the  strong 
detonation  will  be  approximately  a  weak  detonation  followed  by 
a  shock.  The  extreme  situation,  however,  will  occxir  only  for 
extremely  high  reaction  rates. 

Let  us  consider  how  a  detonation  looks  in  the  limiting 
case  when  the  reaction  rate  approaches  zero;   as  far  as  the 
curves  in  the  (  t  ,6,  e  ) -space  are  concer-ad  ""his  is  equivalent 
to  assuming  that  X  and  yu,  approach  zero  since  ^S  A)-  ^^'^f'-'^/Po 
are  the  essential  dlmensionless  parameters.  The  field  vectors 
in  this  case  lie  almost  in  the  planes  £.  =  const,  except  on  the 
lines  cA  and  (@,  As  was  stated  earlier,  the  projection  of  the 
vector  field  in  the  planes  e  =  const, has  a  saddle  singularity  at 
the  points  Bg  ,  Hence  there  are  two  lines  jO+  and  J^^   depending 
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..  For  5=1  they  coincide  with  dO, 
andjOr.   The  curve  <D^      leads  to  larger  values  of  9,  while  ,0" 
leads  to  smaller  values  of  0.   In  the  limit  S  -*»  0,  or,  in  dimen- 
sionless  f orm,  yixS  /p  — >  0  and  XS  /p  -«>■  0,  the  solution  cvLPves 
oP  leaving  the  point  B-,  consist  of  sections  of  the  curve  <S>  followed 
by  the  curves  oO+  or  i5"  until  these  meet  the  plane  r  =  0  or 
the  line  i/r.  Prom  this  remark  it  is  clear  that  the  ultimate 
points  U|3  ~  U/3  (0)  on  the  plane  6  =  0  consist  of  the  cxirve  *(?''"  on 
£  =  0  from  the  point  H  on  the  plane  r  =  0  to  the  point  B  ,  then 
of  the  curve oD"  up  to  the  point  A  .   It  is  clear  from  this  des- 
cription that  the  curve «0*(0)  consists  of  the  curve  &  from  B, 

to  B  followed  by.^"  on  £.=  0  up  to  A^,   In  the  limit  case 

o  *^   o  o   — — ^.^^^— — — — 

/ts/p  -*  0,  Ks/p    -*■  0,  therefore,  the  strong  detonation  consists 

of  a  shock  followed  by  a  reaction  process,  as  expected,-^ 

Finally  we  are  going  to  show  that  weak  deflagrations  exist 

unless  the  reaction  rate  exceeds  a  certain  bound.  To  this  end 

we  should  find  out  whether  or  not  any  of  the  curves  ^   starting 

at  the  point  B^  can  end  up  at  the  point  B  .  The  situation  differs 

from  that  for  detonations  in  that  it  must  now  be  assumed  that 

©  <  0  at  B  ,  (see  p,52  and  Fig,  6),  As  a  consequence  the  line 

'[/(S)  of  ultimate  points  Uo  (S)  is  modified  for  small  values  of  S, 

For  large  values  of  S,  the  situation  is  as  before.  For  S  =  oo  , 

the  initial  part  of  the  ultimate  line1((S)  leads  from  the  point 

H|,  on  'r=  0  straight  over  to  the  i)rojection  of-  H-,  on  &=  0.  For 

large  values  of  S,  therefore,  the  ultimate  lineK(S)  also  begins  at 

H,  and  leads  on  the  plane  T  =  0  over  to  the  plane  e  =  0,  Prom  there 

on  U  leads  on  e  =  0  to  the  point  3^  if  ©  >  0  at  A, ,  or  to  a  point 

U^^  if  ©  <  ©  at  A,  . 
s     1 

For  small  values  of  S,  however,  the  situation  is  quite 
different.  We  can  no  longer  assert  that  the  ultimate  line  leads 
over  from  H,  on  the  plane  T  =  0  to  the  plane  6=0.  As  a  matter 
of  fact,  for  small  values  of  S  the  ultimate  line  stops  on  r  =  0 

1)  If  this  pictixre  of  a  detonation  is  accepted,  then  v,  Neumann's 
result  that  no  weak  detonations  exist  is  implied  by  the  fact  that 
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at  a  point  with  6>  0  and  Jvmps  dls continuously  over  to  the 
line  J^i     the  reason  being  that  a  curve  J5   exists  which  enters 
the  line  (B  below  0  =  ©  ;  (this  possibility  was  excluded  for 
detonations, (see  pp,36,37).  We  first  consider  the  case  S  =0,  In 
that  case  the  curves  J^    consist  of  sections  of  the  c\ipve  (3 
followed  by  sections  of  the  ciarves  £)^   or  SD'   (see  Pig.  4), 
Since  the  curve  (0  now  meets  the  plane  ©  =  ©g  before  &  =  0  is 
reached  it  is  clear  that  none  of  the  curves  J3   will  reach  the 
plane  d-  =  0, 

The  ultimate  line  7^(0)  in  this  case  consists  of  the 
intersections  of  the  lines  JOf  with  r  =  0  from  H^  on  up  to  a 
value  £  =  £g  for  which  d3  intersects  0  =  ©g»  from  there  on 
of  the  part  of  the  curve  t/^with  £  >  £g.  The  ultimate  line 
thus  suffers  a  discontinuity.  Evidently  we  have  e  tf  e^   on  the 
line  7^(0).  This  confirms  the  statement  that  If (0)  does  not  reach 
the  plane  £.=  0, 

It  is  now  clear  that  for  sufficiently  small  values  of  S 
the  ultimate  line  2((S)  will  also  not  reach  the  plane  £=0.  As 
for  S  =  0,  the  ultimate  line  will  consist  of  one  part  on  the 
plane  t  =  0  and  another  part  which  is  a  section  of  the  curve  • 
The  values  of  the  parameter  /3    referring  to  these  two  parts 
are  separated  by  a  value  |S  (S)  such  that  the  solution  curve 
dCj§^  enters  the  point  Gg  for  a  certain  value  e  (S)  of  e  and 
then  enters  the  point  Bg  •   (We  recall  that  Gg  is  the  point  at 
which  the  solution  curve  l^J  through  B^intersects  the  plane  0=0g, 
see  p,35j  For  sujfficiently  large  values  of  S  no  such  value  c 
of  e  exists  since  then  the  ultimate  line  remains  above  ©  =  ©g 
until  it  ends  at  Tf^   or  U^^.  Therefore-  "rhsre  is  a  largest  value 
S^   of  S  for  Tftiich  a  value  e  exists.  It  is  then  seen  that  this 
value  is  e  =  0,  Consequently,  there  is  a  value  /3°  of  /3  such 
that  the  curve  J^°   =  <Oqo  for  S  =  S°  meets  the  plane  £  =  0  at 
the  point  G  and  ends  up  at  the  point  B^.  This  ctirve  J0°   then 
represents  a  weak  deflagration.  Thus  we  have  shown  that  for 
every  value  of  the  flxix  m  and  for  an  appropriate  reaction  rate 
a  weak  deflagration  exists. 
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The  argument  presented  holds  just  as  well  if  the 
points  A,  and  B,  coalesce  so  that  the  flow  corresponding  to 
this  point  is  sonic.  Thus  we  see  that  for  an  appropriate 
reaction  rate  Chapman- Joup:uet  deflagrations  are  possible. 

In  the  limiting  case  where  m  and  S   approach  zero  in 
such  a  way  that  s  =  S^^^  approaches  a  finite  value,  a  defla- 
gration is  always  possible  for  an  appropriate  value  of  S  /m  , 
as  was  shown  in  Section  3,  In  this  limiting  case  dS  =  0  and 
de  =0  for  the  vector  field  except  on  the  plane  9  =  Pt,  on 
which  the  vector  field  is  given  by  (2),  ("3);   (see  p.25).  The 
points  B  and  B-,  lie  on  this-  plane  0  =  p  r  ,  The  line  <£f  ^   is 
the  line  ©  =  ©,  on  e  =  1  with  decreasing  r  .  The  line  ^^  is 
the  line  0  =  Pr  on  e  =  0  with  increasing  r  ,  The  point  G  is 
then  the  intersection  of  9  =  Pr  with  ©  =  0„,  As  is  easily 
seen  there  is  only  one  line  oU  on  the  plane  ©  =  Prand  this  line 
meets  the  point  G^  only  for  a  special  value  s  of  s  =  s/m  , 

Suppose  we  let  m  increase  holding  the  initial  state 
(rQ,©^)  fixed.  Let  S(m)  be  the  value  of  S  for  which  a  deflagra- 
tion exists.  We  have  not  proved  that  S -depends  continuously 
on  m,  but  we  can  be  sure  that  a  continuous  cvirve  of  points  (m,S) 
in  a  (m,S) -plane  exists  representing  pairs  of  values  of  m  and  S 
for  Ti^ich  deflagrations  are  possible.  Eventually  the  fltix  m 
will  reach  a  value  m^  for  which  the  points  A,  and  B,  coalesce 
and  a  Chapman- Jouguet  situation  arises.  Let  S^  be  the  correspond- 
ing value  of  the  reaction  rate  S,  Then  we  can  at  least  say  that 
for  every  value  of  the  reaction  rate  S  5  S«  deflagrations  are 
possible  with  an  appropriate  value  of  the  flux  m. 
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APPENDIX  I 


On  The  Aasmig)tlon  About  the  Discontinuous  Character 
Of  The  Reaction  Process 

The  considerations  of  this  report  rest  on  the  "basic" 
assumption  (see  p.  3)  that  the  reaction  process  may  he  con- 
sidered approximately  a  sharp  discontinuity;  more  specifically, 
that,  firstly,  the  rates  of  change  of  the  pertinent  quantities 
in  the  field  of  flow  outside  of  the  reaction  zone  are  negligibly 
small  when  conpared  to  the  rates  of  change  of  the  same  quantities 
inside  the  reaction  zone  and,  secondly,  that  the  rate  of  change 
of  the  width  of  the  reaction  zone  is  small  when  con5)ar6d  with 
the  average  speed  with  which  the  gases  cross  the  reaction  zone. 
It  was  further  assumed  (p.l9  )  that  the  flow,  when  observed  from 
a  frame  moving  with  the  instantaneous  velocity  of  the  reaction 
front,  is  steady  in  the  neighborhood  of  the  reaction  front  at 
the  time  considered. 

We  first  want  to  show  that  this  latter  assunption  is  con- 
sistent with  the  a3sun5)tion  of  the  discontinuous  character  of 
the  reaction.  Suppose  we  define  the  velocity  x  of  a  border  of 
the  reaction  zone  by  u^  +  xu^  =  0,  and  suppose  we  choose  the 
velocity  of  our  frame  such  that  x  =  0  at  one  point  inside  the  re- 
action zone.   Then  the  assximption  that  the  rate  of  change  of  the 
width  of  the  reaction  zone  is  small  compared  with  the  average  gas 
velocity  In  the  zone  can  then  be  fornailated  as  |i|«  u  everywhere 
in  the  zone.   Consequently,  |u^|«|uu^!  -  rywhere  in  the  zone. 
The  term  u. ,  can,  therefore,  be  omitted  from  the  differential  equa- 
tion which  expresses  the  fact  that  the  acceleration  u.  +  uu  equals 
the  total  applied  force  per  unit  mass.  For  similar  reasons  one  can 
omit  the  terms  r. ,  0. ,  and  £^  from  the  differential  equations 
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expressing  the  balance  of  mass  flow,  energy,  and  chemical  re- 
action. In  other  words,  to  the  degree  of  accuracy  Implied  by 
our  basic  assumption,  the  equations  characterizing  non-steady 
flow  reduce  to  the  equations  (1)',  (2)',  (3)   In  addition  to 
m  =  const.   Thus  our  assumption  of  "local  steadiness"  Is  in 
agreement  with  o\ir  "basic"  assumption. 

Secondly  we  want  to  mention  that  frequently  reaction 
flow  processes  occur  In  which  our  basic  assumption  is  not 
satisfied.  Detonations,  consisting  of  a  chemical  reaction 
process  initiated  by  a  shock,  are  frequently  followed  by  rare- 
faction waves.   It  is  clear  that  this  rarefaction  wave 
interferes  with  the  reaction  process,  but  the  interference  can 
be  ignored  if  our  basic  assumption  is  satisfied  and  the  changes 
which  the  rarefaction  wave  produces  in  a  section  of  the  width 
of  the  reaction  zone  are  insignificant.  If,  however,  the  re- 
action zone  is  so  wide  that  this  interference  can  no  longer  be 
ignored,  then  pressure  and  teii5)erature  in  the  reaction  zone  are 
diminished  and  the  strength  of  the  Initiating  shock  is  reduced. 
In  particular,  the  speed  of  the  detonation  wave  is  then  leas 
than  that  calculated  without  interference  and  could  thus  be 
less  than  that  of  a  Chapman- Jouguet  detonation.  If  the  inter- 
action is  strong  the  detonation  may  eventually  cease. 

The  occurence  of  a  lower  detonation  limit  may  be  ex- 
plained in  this  way,  since  the  reaction  rate  is  low,  and  hence 
the  reaction  zone  is  wide,  if  the  concentration  of  the  com- 
bustible component  in  the  unburnt  explosive  mixture  is  low. 
(See  Wendlandt  [13] ). 

Spherical  detonation  waves  maintaining  the  Chapman- 
Jouguet  condition  are  possible  inasmuch  as  the  detonation  pro- 
cess can  be  considered  a  discontinuity,  as  Taylor  has  shown 
[14].  The  rarefaction  wave  following  the  detonation  front  be- 
gins, however,  with  an  infinite  rate  of  change  of  the  signifi- 
cant quantities,  (see  [14]  )•      Consequently,  this  rarefaction 
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wave  always  Interferes  with,  the  detonation  wave  noticeably  and. 
diminishes  its  strength  and  speed.  It  can  thus  be  understood 
why  the  speed  of  spherical  detonation  waves  Is  less  than  re- 
quired by  the  Chapman- Jouguet  hypothesis.   Calculations  of  such 
spherical  detonation  waves  with  sub-Chapman-Jouguet  speed  were 
carried  out  by  Eyrlng  and  his  collaborators  [6],  [V]  . 

The  gradual  building  up  of  a  detonation  wave  ignited  at 
the  closed  end  of  a  tube  may  also  be  explained  as  due  to  the 
interference  of  a  rarefaction  wave  with  the  reaction  process. 
The  condition  that  the  flow  velocity  vanish  at  the  closed  end 
requires  that  a  sirple  rarefaction  wave  follow  the  detonation 
wave.   This  rarefaction  wave,  when  it  begins,  involves  an  in- 
finite drop  of  pressure,  temperature  and  velocity;  thus  a 
noticeable  Interference  is  clearly  Indicated.   (For  numerical 
calculations  of  the  gradual  building  up  of  a  detonation  see  [7]  )• 

As  regards  deflagration  processes  it  would  seem  probable 
that  conibustion  processes  occur  in  which  the  width  of  the  re- 
action zone  widens  noticeably.   Otherwise  It  would  seem  im- 
possible to  explain  how  a  flame  could  ever  overtake  a  shock  front 
preceding  it.  Such  combustion  processes  would  then  be  analogous 
to  rarefaction  waves  and  not  to  shock  discontinuities.   No 
theoretical  treatment  of  such  processes  seems  to  exist. 

These  remarks  are  Intended  to  show  that  the  frequently 
observed  deviations  from  the  predictions  of  the  discontinuity 
thoory  are  not  due  to  the  unsteadiness  of  the  process  as  such 
but  rather  to  the  occurrence  of  a  relatively  wide  reaction  zone 
which  permits  the  Interference  of  the  outside  flow  with  the 
reaction  process. 
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APPENDIX  II 


It  was  shown  in  the  text  that  for  an  "excessively" 
high  reaction  rate  S  =  S(p,©)  a  weak  detonation  occurs  In- 
stead of  the  Chapman- Jouguet  detonation  (see  p.  38  ).   It 
Is  of  interest  to  know  how  high  a  reaction  rate  must  be  in 
order  to  be  "excessive"  in  this  sense.  The  excessive  cases  are 
separated  from  the  regular  ones  by  "maximal"  cases  In  which  a 
Chapman- Jouguet  detonation  is  just  possible,  while  such  a  detona- 
tion is  impossible  for  a  reaction  rate  higher  than  a  maximal  one. 
We  shall  present  some  such  maximal  reaction  rates  nximerlcally. 
Prom  these  results  it  will  appear  that  for  maximal  reaction 
rates  the  transition  zone  becomes  extremely  small,  of  the  order 
of  magnitude  of  one  mean  free  path,  if  viscosity  and  heat  conduc- 

Under 

these  circiimstances  the  notions  of  viscosity  and  heat  conduction 
in  the  transition  zone  become  meaningless.  If,  however,  viscosity 
and  heat  conduction  are  ten  times  as  large  as  for  atmospheric  air 

For  ex- 
>  excessive  if  it  vanishes  up  to  900°I 
equals  16,10   sec"  or  more  for  higjher  temperatxires ,  The  width 
of  the  transition  zone  wovild  then  be  about  ten  times  as  long  as 
a  mean  free  path. 

We  assiome  that  unburnt  and  burnt  gas  are  polytroplc  with 
exponents  ^  =  1,4  and  i  =  1,2,  For  the  mixture  consisting  of 
the  fraction  e  of  burnt  gas  and  1- t    of  unburnt  gas  we  define  V 

by 

1   «   -  '-  £  +   e 
V-l~   Vo-1   Vi-1   • 

Then  we  have  for  the  energy  per  vmlt  mass  of  the  mixture  the  ex- 
pression 

E(^J(e)=— ^  +  (1  -£)P   , 
V-  1 
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in  which  the  liberated  energy  per  unit  mass  f  is  assumed  to 

be  _ 

r  =  32.5  ©Q 

which  corresponds  to  a  value  of  665  cal/gm  if  the  initial  tempera- 
ture is  M  e  /R  =  300°K  and  the  molecular  weight  of  the  xmburnt 
o  o'  o 

_,  =  29. 

About  the  heat  conductivity   X  and  viscosity^  we  have  made 
the   assumption  ^ 

/^~  y-1 

which  was  proposed  by  Becker  [4].  We  shall  in  particular  consider 

as  reference  value  forjx  the  value 

-4   / 
lA.     =  2«44»10  gm/cm  sec  . 

with  ^i=     V  =  1,4,  the  corresponding  value  of  X  is  then 

\  =  8»54»10"  gm/cm  sec. 


We  recall  that  the  customary  coefficients  of  viscosity  and  heat 

4 
3 


conductivity  are  ^lx,    and  RX  in  our  notation;   (see  footnote  3  on 


p.l2). 

The  reaction  rate  S  was  assumed  to  be  zero  up  to  a  safety 
temperature  ©  =  ©„;  for  0  >  Og  we  have  asstimed 

iU  S  =  const ,  ; 
thus,  if  u  is  independent  of  temperature  and  pressure  the  same 
is  then  assumed  of  the  reaction  rate  for  Q  >  Oq.  This  assumption 
is  rather  imrealistic  because  both  the  reaction  rate  and  the  vis- 
cosity will  increase  with  the  pressure;  but  this  asstmiption 
should  be  sufficient  to  given  information  about  the  order  of  mag- 
nitude of  maximal  reaction  rates.   In  Pig.  6  we  have  plotted 
such  maximal  reaction  rates,  or  rather  the  values  -^  S,  in  which 
^  is  the  reference  viscosity  given  above.  For  a  s^ety  tempera- 
ture of  900°K  or  0  s  3©  ,  for  example,  we  find  as  maximal  value 


51 


i.' 

171  sec~ 

lo" 

/ 

/ 

/ 

lo" 

/ 

li" 

^  ^ 

/ 

/ 
/ 

Pig.   6 
Maximal  reaction  rates  S   as  fimctions  of  the  temperature  T, 


(If  the  reaction  rate  exceeas   a  maximal  one,    a  Chapman- 
Jouguet  detonation  is  not  possible, /^y^    is  the  ratio  of 
actual  viscosity  to  that   for   atmospheric   air   at  room 
temperature.) 

q    /i  -% 

Of  the  reaction  rate  S  =  16.10     ~  sec"     for  ©  >  ©„. 

jtoi 

unburnt  gas  at  rest,  solely  determined  by  the  Chapman-Jouguet 

condition,  equals  v  =  15*7  msec"  ,  The  width  of  the  reaction 

zone  is  roughly  given  by 


S  =  v^s- 


-1 


In  the  case  ©  =  5©  ,  S  =  16»10  <-^   sec""^,  we  therefore  have 

S  =  98.10"'''  -^   cm  , 

a  length  which  for  j^=     M-  would  be  about  equal  to  one  free  mean 
path  in  air. 
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To  determine  maximal  reaction  rates  we  observe  that 
the  Chapman-Jotiguet  condition  cannot  be  satisfied  if  it  would 
lead  to  a  situation  called  "abnormal"  in  Sec,  5, p. 34.  In  ab- 
normal situations,  the  curve  ^  starting  at  the  point  v/^,  ends 
up  on  the  plane  d  =  0  with  a  value  ©  >  ©  . 

A  maximal  situation,  in  which  the  Chapman- Jouguet  condi- 
tion can  just  be  satisfied,  therefore  corresponds  to  ^at  was 
called  an  "extreme"  situation,  in  which  the  cvtrve  ^  ends  up  on 
5=0  just  with  ®  =  ®s»  To  obtain  a  maximal  reaction  rate  we 
then  proceed  as  follows.  To  a  given  initial  state  (t'o'Pq^o^ 
corresponding  to  the  point  ^/V  ,  we  determine  the  end  state  (''vP^ ,  ©,) 
from  the  Chapman- Jouguet  condition.  The  flux  m  is  then  also  de- 
termined. We  now  assume  any  value  for  the  reaction  rate  S,  or 
rather  for  uS  and  determine  the  curve  ^.  The  value  of  ©  with 
which  ^  ends  upon  £.  =  0  is  then  taken  as  safety  temperatiire  ©  • 
We  finally  select  those  values  for  the  reaction  rate  S  for  which 
the  safety  temper at\ire  turns  out  to  be  between  600*^K  and  900°K, 

The  curve  ^  was  characterized  as  the  graph  of  that  solu- 
tion of  the  differential  equation  which  leaves  the  point  i/V,  with 
decreasing  x.  All  other  solution  curves  entering  J^^   lie  on  £=  1, 
Consequently,  ^  can  also  be  characterized  as  the  only  solution 
curve  leaving  the  point  v/V  without  de  =0,  Hence  &   can  be  intro- 
duced as  parameter.  The  differential  equations  then  become,  after 

2 
introducing  t  =  m  T  , 

(1-6) /AS  11  =  t(p+t-  P-^-t^), 
(l-£)VyiiS  ^=t[pt  -^i-t^  +  (V-l)tt^  +  (?-l)(l-£)fPit3^-Plt^-^t2], 

The  desired  solution  t  =  t(£),  p  =  p(e)  is  then  the  one  that  assumes 
the  values  (t,,p,)  for  £.=  1  and  permits  expansion  with  respect  to 
powers  of  (l-£).  It  is  easily  obtained  from  this  power  series 
for  small  values  of  (l-£)  and  by  finite  differences  for  larger 
values  of  1  -  £  up  to  £  =  0, 
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